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ABSTRACT 



A computer code was developed for the analysis of air- 
cooled heat exchangers and was coupled with a constrained 
function minimization code to produce an automated air-cooled 
heat exchanger design and optimization program with many 
new capabilities. 

A general iteration free approximation method was used 
for the analysis which calculates the mean overall heat 
transfer coefficient and the overall pressure drop for many 
flow arrangements, taking into account the variation of the 
heat transfer coefficients and the pressure drop with tempera- 
ture and/or length of flow path. 

The code is not limited to surfaces found in the litera- 
ture, but will accommodate any triangular pitch bank of finned 
tubes in multiple-pass configurations. 
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NOMENCLATURE 



English Letter Symbols 

2 

A - total heat transfer area, in 

0 ^ - specific heat, BTU/lbm-°F 

C - heat capacity rate, BTU/hr-°F = m c^ 

D - diameter, in. 

f - friction factor 

F - LMTD correction factor 

2 

g c - acceleration of gravity, 32.2 ft/sec 

2 

H - corrected heat transfer coefficient, BTU/hr-ft -°F 
h - bank height, in. 

J - Colburn factor 

k - thermal conductivity, BTU/hr-ft- °F 

L - length, in. 

l - fin height, in. 

m - mass flow rate, lbm/hr 

m - / 2H /k,t , ft -1 

o f 

n - number of (used with appropriate subscript) 

N - number of tubes 

p - pressure, psi 

P - pitch, in. 

Pr - Prandtl number 

£) - heat transfer rate, Btu/hr 

r - radius, in. 

2 

R - heat transfer resistance, hr-ft -°F/BTU 
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Re - Reynolds number 

s - distance between adjacent fins, in. 

S - fin spacing center-to-cent er , in. 
t - fin thickness, in. 

T - true temperature, deg. 

2 

U - overall heat transfer coefficient, BTU/hr-ft -°F 
w - bank width, in. 

Greek Letter Symbols 

AT - temperature difference, deg. 

n - surface efficiency 

cp - fin efficiency 

y - viscosity, lbm/ft-hr 

3 

p - density, lbm/ft 

- temperature correction 

Subscripts 



a 


— 


axr 


c 


- 


cold 


f 


- 


fin 


ff 


- 


free face 


h 


- 


hot 


i 


- 


inside 


j 


- 


reference number 


l 


- 


limiting 


L 


- 


longitudinal 


m 


- 


mean 


o 


_ 


outside 
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p - pass 

r - rows 

t - transverse 

T - total 

w - water 

x - cross-sectional 

1 - entering 

2 - leaving 

®> - ambient 

Superscripts 

uncorrected 
* - estimated 

i - initial 



(flow) 
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I . INTRODUCTION 



A. BACKGROUND 

The cooling of fluids by passing ambient air over extended 
tube surfaces is a relatively recent development in heat 
exchangers. Its application has come about cautiously, due 
to the usual reluctance to change from well established and 
well documented methods, i.e. the shell- and- tube heat exchanger. 

However, concern for the environment and economic pres- 
sures have necessitated the use of air as a coolant. Smith 
[1] has listed some typical advantages of direct cooling with 
air as compared to cooling with water in a shell-and-tube 
exchanger : 

a. Eliminates the problem of temperature rise in, and 
pollution of, water resources. 

b. Enables plant location to be independent of a water 
supply . 

c. Eliminates the necessity of much coolant piping. 

d. Reduces heat exchanger maintenance costs by eliminating 
the need of descaling water-side surfaces. The 
mechanical drives will operate in a noncorrosive 
atmosphere . 

e. Eliminates water treatment. 

f. Limiting coolant temperatures is unnecessary. 

g. Enables installation of exchangers at elevations above 
other operating equipment at no penalty, thus reducing 
ground area requirements. 

With air colling becoming more and more competitive with 
water; even when water supplies are plentiful, an automated 
air-cooled heat exchanger design package (which could be used 
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for trade-off studies, first cut analysis, and conceptual 
design), would be of great use. 

The design of an air-cooled cross flow heat exchanger is 
a complex task requiring the examination and optimization of 
a wide variety of heat transfer surfaces. Studies have shown 
that a poor choice of either the heat transfer surfaces or 
design parameters can more than double the costs charageable 
to a heat exchanger [2], 

For the optimized design of heat exchangers with the 
computer; reliable, but fast, calculation methods for the 
mean overall heat transfer coefficient and the overall 
pressure drops are needed for the following reasons; 

a. Conventional simple methods using mean values of 
temperatures as reference temperatures can lead to 
undesirable errors [3] . 

b. Numerical stepwise integrations are prohibitively 
time consuming . 

B . REVIEW 

A number of heat exchanger design methods have been proposed 
to determine the optimum heat exchanger design. Bergles, et 
al. [4], performed an evaluation of different objective func- 
tions for compact heat exchangers with different heat trans- 
fer surfaces, but the same specifications. The method did not 
include any actual optimization techniques, but results did 
show that a great improvement in heat exchanger performance 
can be made by proper selection of design parameters. 
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The method of Fax and Mills [5] , used Lagrange multi- 
pliers to optimize a heat exchanger design under specified 
constraints. This technique required that the objective 
function and constraints be expressed explicitly and be differ- 
entiable throughout the range of interest. The total num- 
ber of constraints had to be less than the total number of 
variables, and all constraints had to be equality constraints. 
The method was obviously restricted to a very limited number 
of problems. 

For a unit section of an air-cooled exchanger of standard 
length and width equipped with specified fans, both Schoonman 
[6] , and Joyce [7] , used factorial searches of fin spacing, 
number of rows and air rate to maximize the ratio of heat 
transfer to cost. Because the number of sections must be 
rounded to an integer, this method only appears suitable for 
large exchangers . 

Nakayama [3], used a similar approach requiring the plotting 
of heat transfer coefficients and pressure drops which pre- 
sumably could be programmed into a factorial search. 

Kern [8] derived an analytic expression for the annual 
cost as a function of air rate and number of rows. (This was 
done by using a constant tubeside heat transfer coefficient 
and an arithmetic mean temperature difference.) The optimum 
was found by setting partial derivatives equal to zero. 
Alternatively, if the airside is assumed to control, the 
optimum allocation of the components of total cost might be 



14 



found through geometric programming, as illustrated by Auriel 
and Wilde [9]. Oshwald and Kochenberger [10], also presented 
a geometric programming method for heat exchanger optimization 
and used it to select heat exchanger fluids considering power 
requirements, cost, tube diameters, velocities, temperature, 
and other physical properties. 

After discussing factorial, univariate, and random search 
methods for an optimum design of a shell-and-tube exchanger, 
Briggs and Evans [11, 121, discuss a "logical search method"; 
or what Peters and Nicole [13], call "heuristics". With this 
method, an engineer makes use of selecting design variables 
close to the optimum to obtain an optimum design. The heuris- 
tic method is less scientific, and is useful only when computer 
time and storage are at a premium. Due to the large number 
of discrete variables encountered in air-cooled heat exchanger 
design, Peters and Nicole [13] , chose to base their cost- 
optimizing design programs on heuristic algorithms (starting 
close to optimum) , specific to the equipment under considera- 
tion. 

Mott, et al. [2], discuss a computerized procedure for 
designing a minimum cost heat exchanger. The method minimizes 
a cost index expressed as a function of fluid pumping power. 

The algorithm imposed no constraints. 

To this point, no mention has been made of applying the 
concepts and techniques of nonlinear programming to optimizing 
the design of heat exchangers. However, Palen, et al . [14] , 

in 1974 proposed using the Complex Method [15] , for the heat 
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exchanger optimization problem. They found a minimum cost 
shell-and-tube exchanger by varying six geometrical parameters. 
The Complex Method requires several feasible starting designs 
before optimization can be performed. 

Johnson, et al . [16], coupled an existing shell-and-tube 

condenser design code with a constrained function minimization 
code to produce an automated marine condenser design program 
of vastly different complexity. 

The most complete work to date has been accomplished by 
Afimiwala [17] . He has applied various nonlinear programming 
methods of optimization to the heat exchanger design problem; 
including an experimental interactive graphical approach and 
exterior penalty function techniques. The gradient based 
search methods of Davidson - Fletcher - Powell and conjugate 
gradient were used for the resulting unconstrained minimiza- 
tions. The exterior penalty method is extremely useful, since 
an initial solution satisfying the constraints is not required. 
The gradient based search methods are efficient when considering 
computer time. 

Finally, Fontein and Wassink [18] utilized the complex 
method of Nelder and Mead [19] , and a steepest descent method 
[20] , for optimizing a shell-and-tube exchanger. 

It can be seen that although there are many methods that 
have been presented for heat exchanger optimization, each of 
the methods has its own limitations; none is completely general. 
Of all the design procedures cited above (those of which are 
applicable to cross flow air-cooled heat exchangers) , all are 
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limited to the 120 individual surfaces found in the open 
literature [21] for the calculation of the air-side heat 
transfer coefficient and friction factor. Therefore, the 
designer is faced with choosing an optimum surface from a 
number of individual optimal designs calculated from one of 
the above methods. In addition, the above methods treat the 
overall heat transfer coefficient as a constant, or they become 
involved with time-consuming numerical stepwise integrations 
in an attempt to account for the varying heat transfer coeffi- 
cients . 

This paper tries to bridge this gap by presenting an 
optimization routine that: selects an optimal surface, takes 

into account the varying heat transfer coefficients and fric- 
tion factors across the exchanger, performs each analysis in 
an iterative-free manner, and may start with an infeasible 
design . 

C. METHODOLOGY 

With the Control Program for Engineering Synthesis and 
Con strained Function Min imization (COPES/CONMIN) optimizing 
scheme, a nonlinear optimization program is available that is 
capable of optimizing a wide class of engineering problems 
[22,23]. Therefore, for the heat exchanger design problem, 
it was necessary to develop a subroutine, which given a start- 
ing design, would analyze an air-cooled heat exchanger, and 
which would be suitable for coupling with the optimizer, 
COPES/CONMIN. 
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Figure 1 illustrates the procedure by which the heat 
exchanger was analyzed. Initial input consisted of a com- 
plete listing of all design parameters, whether known or 
estimated. Those that were estimates, i.e. unknowns, were 
later passed to the optimizer as design variables. 

The analysis scheme and the optimizer will be discussed 
in much greater detail in the following chapters. 

D. OBJECTIVES 

The objectives of this thesis are two-fold. 

The first objective was to develop a computer code, 
hereafter referred to as ANALIZ, which would analyze an 
air-cooled heat exchanger given any initial design. The 
analysis scheme was to: be iterative-free, take into account 

the variation of the heat transfer coefficients and the pres- 
sure drop with temperature and/or length of flow path, and 
finally, be written in such a manner that it could be coupled 
with an existing optimizer COPES/CONMIN. 

The second objective was to actually couple ANALIZ with 
a numerical optimization program. This would produce a 
detailed design program which would have the capability 
to determine an optimum surface, while optimizing the objec- 
tive function such as size, weight, cost, etc. 
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II. NUMERICAL OPTIMIZATION 



A. BACKGROUND 

Almost all design problems require either the maximiza- 
tion or the minimization of some parameter or function. 

This parameter shall be called the design's objective function 
[24] . For example, the problem may call for a heat exchanger 
with a minimum volume. The expression for volume would be 
the design objective function. For the design to be accepta- 
ble, it must satisfy certain design constraints. For exam- 
ple, an air heater must be designed so that it will fit into 
a given space. Therefore, the engineer must set design con- 
straints on the maximum size of the exchanger. 

If the objective function could be easily formulated 
analytically, the maxima or minima could be found by using 
the methods of differential calculus. However, the limita- 
tions of this method are obvious . 

Another numerical method that would be satisfactory for 
small scale problems would be an iterative solution technique. 
A computer program could be written containing a series of 
nested iteration loops that would vary the design parameters 
and solve the problem for a variety of values for each of the 
parameters. For other than small, easily formulated problems, 
the cost in central processor (CPU), time would be prohibitive. 

Over the last twenty years, many numerical optimization 
techniques have been developed specifically for computer 
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utilization. These techniques usually do not require a 
specific algebraic equation, but rather any computer algorithm 
to which design variables can be input and from which the 
objective function and design constraint values can be 
determined is acceptable. For this reason, nonlinear pro- 
gramming methods were chosen for the air-cooled heat exchanger 
design. Some of these techniques were summarized by Shah et 
al . [25] in figure 2. 

1. One-Dimensional Search Methods . Two of the most common 
of these search methods are the golden section [2.6] 

and quadratic interpolation [27], The former iso- 
lates the minimum in regions of successively decreasing 
size, the latter performs a series of iterations 
approximating the objective function as a quadratic . 

2. Multidimensional Unconstrained Search Methods . These 
unconstrained searches can be performed by a sequence 
of one-dimensional minimizations in the proper direc- 
tions . 

3 . Multidimensional Constrained Search Methods . A common 
method here for enforcing the constraints in an optimi- 
zation scheme are based on the sequential penalty func- 
tion method. These techniques convert the constrained 
optimization problem into a sequence of unconstrained 
problems. This is accomplished by applying either an 
exterior or interior penalty to the objective function. 
The Complex Method [15] , locates the optimum based on 
an intuitive approach in a n-dimensional space defined 
by the independent design variables. The method of 
feasible directions is used primarily for inequality 
constraints . 

An optimization program based upon the Augmented Lagrangian 
Multiplier Method and the method of feasible directions was 
chosen for this research project. 



B. CONSTRAINED FUNCTION MINIMIZATION (CONMIN) 

Vanderplaats [22], developed an optimization program, 
CONMIN, based on the method of feasible directions which is 
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capable of optimizing a wide variety of engineering problems. 
CONMIN is a FORTRAN program, in subprogram form, that opti- 
mizes a function subject to a set of inequality constraints. 

The following definitions will be useful in the following 
discussion : 

1. Design Variables - those parameters which the optimiza- 
tion program can change in order to improve the design. 

2. Design Constraints - those parameters which must not 
exceed given bounds for the design to be acceptable. 

3. Objective Function - the parameter which is going to 
be minimized or maximized. 

The general nonlinear inequality constrained optimization 
problem can be written mathematically as follows [28]: 

( 1 ) 

0 j = 1 , NCON (2) 

< X U k= 1 , NDV (3) 

— l 

(4) 



The vector X is the vector of design variables, with NDV 

equal to the number of design variables. The objective function. 



where 



Minimize F (X) 
Subject to: 

g j (x) i 



XV < X. 

1 — 1 



x. 



X = 



X, 



X 



NDV 
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F (X) , given by eq . (1), as well as the constraint functions 

given by eq. (2) , may be linear or nonlinear functions of 
the design variables. They shall also be explicit or implicit 

functions of X, but must have continuous first derivatives. 

£ u. 

NCON is the number of constraints. X. and X. are the lower 

1 i 

and upper bounds or side constraints placed on the design 
variables. Side constraints could be included in eq. (2), 
but are treated separately for efficiency. Equality con- 
straints are not dealt with by CONMIN, but will be treated 
separately by a multiplier method. 

CONMIN requires that an initial vector of design variables, 

X, which may or may not yield a feasible design, be specified. 

The design process continues iteratively as: 

3c q+1 = + a* S* 3 (5) 

where is a vector search direction, a* is a scalar quantity 
which determines the amount of change in X and q is the itera- 
tion number. At iteration q a direction must be found 
which will reduce the objective (usable direction for mini- 
mization) , without violating any constraints (feasible direction) , 
see Figure 3 [30] . Once is determined, eq. (5) becomes a 
one-dimensional search problem in whcih a* must be found 
such that; F(X) is at a minimum (see Figure 4 [29]), a new 
constraint is encountered, or a currently active constraint 
(gj (X) = 0) is encountered again. 

The design problem at iteration q+1 becomes one of 
finding a usable-feasible direction, S^, and a move parameter 
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a*. This process is illustrated geometrically by Johnson 
[30] . Consider a condenser problem with just two design 
variables, and X^, where 



X^ = condenser tube outside diameter, 
X 2 = tube pitch to diameter ratio. 



Let the objective function be condenser volume, VOL (X) . 

Assume that the tube bundle diameter must be greater than a 

given value, BD min ' and that the cooling water pumping power 

must be less than a given value HP . Figure 5 illustrates 

3 max 

the design problem geometrically. 

It should be reiterated here, that while Johnson's exam- 
ple starts with a feasible initial design, (a) , this is not 

a requirement and CONMIN is capable of optimizing given an 
infeasible initial design. This is obviously a great benefit. 

The optimization process begins by calculating the gradient 
of the objective function by using finite difference. Each 
design variable is perturbed by .01 in a single forward step. 

The gradient of the objective function, VF , shown in 
Figure 6 is simply the vector of the first partial deriva- 
tives with respect to the design variables; that is: 



VF(X) e v VOL = 



3F 




AVOL 


3X. 




AX, 


1 


• 


1 


3F 




AVOL 


3X 2 




ax 2 



( 6 ) 
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Therefore, because no constraints are active or violated at 
(X) . VF defines the direction of steepest ascent. Because 
it is desired to minimize F, the greatest improvement can be 
made by moving in the negative gradient direction so that 

S = - VF = - VVOL . (7) 



With the value of S now determined, a search is per- 
formed from (X) until the minimum F is found at (X^) on Figure 
6. This is accomplished by taking several values of in 

eq. (5) and interpolating for the a* which will give the mini- 
mum value of F. 

The second design iteration is begun at (X^ by again 
perturbing X to find VF . Instead of moving in the direction 
of steepest descent, a conjugate direction, developed by 
Fletcher and Reeves [31] , is chosen with this method, s’ is 
calculated as follows: 



_ _ lVF(X) q | 2 

S* = - ’F(X) q + lv P (jf,q-l|2 



: q-l 



( 8 ) 



see Figure 7 . 

The Fletcher- Reeves method is used in order to speed 
convergence. With the new conjugate direction, a search is 
performed in this direction until a constraint is encountered. 
This occurs at © of Figure 8 on the pumping power constraint. 

At (X) , with the HP constraint active, not only is VF 
found, but the gradient of the active constraint is also 
computed, again using finite difference. The requirements on 
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the new search direction are now twofold; it must reduce the 
objective function and, at the same time, not violate the 
active constraint. This is solved by using the method of 
feasible directions developed by Zoutendijk [32] , and imple- 
mented by Vanderplaats and Moses [33] . 

The problem of finding the new S can be stated as [34] : 

Maximize 6 subject to the constraints 



VF (X) * 


S + g <_ 0 




(9) 


(X) 


•s + 0 . g < 0 

3 


j = 1 , NAC 


(10) 


s ♦ s 


< 1 




(ID 



where Vg^ (X) = - "VHP and NAC is the number of active constraints 
(in this case NAC = 1) . 

If equation (9) is satisfied and g is positive, the search 
direction will reduce the objective function and is defined 
as a usable. If equation (10) is satisfied and g is positive, 

S is a feasible direction, because in this direction, no con- 
straints will be violated if only a small move is taken. See 
Figure 3. 0^ is defined as the push-off factor for the active 

constraint and causes the design to move away from the con- 
straint. Gj must be greater than or equal to zero in order 
to maintain a feasible design. If the maximum value of g from 
equations (9) through (11) is zero, then there is no direction 
that will both reduce the objective function and also be feasi- 
ble. Therefore, the current design is, at least, a local 
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minimum. In Johnson's example, a usable-feasible direction 

exists and a one-dimensional search leads to (j}) in Figure 8 

where the minimum bundle diameter, BD . , constraint is met. 

min 



From 
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it should be noted here that CONMIN had informa- 



tion regarding the linearity of the BD . constraint and, 

min 

therefore, in (10) , has set = 0 to allow S to follow the 
constraint as shown in Figure 9. The one-dimensional search 
along this constraint is carried out until no further design 
improvement is realized. This occurs at . 

This discussion of CONMIN would not be complete without 
citing the program's limitations. NDV directly affects the 
computational time required to reach the optimum. Since the 
calculation of gradients required for each design variable at 
the beginning of each design iteration is found by a finite 
difference step, which requires a complete pass through the 
analysis portion of the program, there is a subsequent in- 
crease in CPU time as NDV increases. Also, due to the inter- 
action between design variables, as NDV gets larger, convergence 
slows during the optimization process. Vanderplaats [24], 
recommends that for most problems of general interest a prac- 
tical limit of NDV = 20 be imposed. NCON does not present 
the same problem because gradient information is calculated 
simultaneously with VF and then only if the constraint is 
active or violated. 

CONMIN offers no guarantees that a global minimum has been 
reached. Therefore, to lend some assurance, the design is 
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started with several different initial vectors until the 
same optimal design is reached. 

Although CONMIN performs very well with inequality con- 
straints, equality constraints such as: 

h R (X) = 0 

cannot be dealt with directly, but must be treated separately, 
using a different method, which will be discussed in the 
following sections. 

C. CONTROL PROGRAM FOR ENGINEERING SYNTHESIS (COPES) 

Recall that CONMIN was written in subroutine form, 
Vanderplaats [23] , has developed a main program which greatly 
enhances the use of CONMIN. 

For this main program, COPES, the user must supply an 
analysis subroutine titled ANALIZ. The subroutines CONMIN 
and ANALIZ are then used by COPES to optimize the objective 
function subject only to the inequality constraints. 

ANALIZ must be organized into three segments: input, 

analysis and output. Based on the value of a counter, ICALC , 
ANALIZ performs the proper function in sequence. 

The COPES program currently provides four specific capa- 
bilities : 

1. Single Analysis - one cycle through the program, as 
if ANALIZ was executing alone. 

2. Optimization - minimization or maximization of the 
objective function with constraints and side- 
constraints imposed. 
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3. Sensitivity Analysis - used to explore the effect of 
changing one or more design variables on one or more 
functions . 

4. Two Variable Function Space - provides tables of data 
of all specified combination of two design variables. 

However, a recent addition to COPES (still in the develop- 
mental stages) , has put the Augmented Lagrangian Multiplier 
Method (ALMM) at the disposal of the programmer. Because of 
its good rate of convergence and its theoretical properties, 
the ALMM is preferred for equality constrained problems [35] . 

Therefore, COPES will call the various subroutines in order 
to optimize the objective function; subroutine ANALIZ for 
necessary analysis information, CONMIN for an optimum based 
on inequality constraints, and a subroutine (yet to be named), 
utilizing ALMM for an optimum satisfying the equality con- 
straints . 

For detailed explanation of the ALMM, its background and 
mathematical derivation, consult reference [35]; but for now 
consider the equality constrained problem; 

Min f (X) (10a) 

Subject to h k (X) = 0 k = l,NECON (10b) 

where NECON is the number of equality constraints. 

Define the "modified" Lagrangian function as; 

m 

L (X , A ) = f(X) + l A. h. (X) 

i=l 1 K 

where A ^ is the Lagrangian multiplier. The problem can now 
be stated as: 
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Min L (X , A ) 



(10c) 



Subject to h^(X) = 0 (10d) 

Then, according to Lagrange, if a I can be found for which 
X solves the problem stated above, then X is also the solution 
to the original problem, eq. (10a) and (10b) . 

The new problem is solved by the conventional exterior 
penalty function method, because this is believed to be one 
of the most efficient algorithms for the solution of such 
equality constrained problems [35] . 
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II. HEAT EXCHANGER ANALYSIS 



A. INTRODUCTION 

To meet the objectives of this thesis, an analysis program 
for an air-cooled heat exchanger must be coupled with a 
numerical optimization scheme to produce a complete, detailed 
design package. COPES/CONMIN has greatly simplified this task 

This analysis program must be written in subroutine form, 
titled ANALIZ, and organized into three segments: input, 

execution, and output. The analysis subroutine must also: 

1. Take into account the variation of the heat-transfer 
coefficients and differential pressure drop with 
temperature and/or length of flow path. 

2. Be iterative free, if possible. 

3. Be written in such a manner that the optimizer will 
play a role in surface selection. 

With the number of design variables approaching the prac- 
tical limit, the importance of an iterative free analysis 
subroutine cannot be over-emphasized. The reason being, that 
at the beginning of each design iteration in CONMIN, the calcu 
lation of all gradients (each design variable and active con- 
straint) requires a complete pass through ANALIZ. Therefore, 
the computational time required by ANALIZ directly affects the 
time required to reach the optimum. 

B. PROBLEM FORMULATION 

The air-cooled heat exchanger is shown in Figure 10. A 
cross-flow arrangement with both fluids unmixed was chosen. 
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Cool air enters the heat exchanger at temperature T , 

C 1 

pressure p^, and constant specific heat c . The cool air 

p a 

makes one pass through the exchanger as it flows over an isos- 
celes pitched bank of finned tubes. The air is heated by 

water, in single phase, at an entering temperature of T, , 

n l 

and constant specific heat c 

p w 

The analysis of the air-cooled , cross-flow heat exchanger 
centers about the first law of thermodynamics and on the heat 
transfer equation. These equations as they apply to the 
exchanger of Figure 10 are as follows: 



Q 0 = m c (T - T ) 
3 a p a C 2 C 1 



( 12 ) 



Q A = m e (T. - T. ) 

4 w Pw n l h 2 



(13) 



Q c = U A AT 
5 mm 



(14) 



where m , T and m , T, are the fluid mass flow rates and 

a c 2 w h 2 

exit temperatures of air and water respectively. is the 

true mean overall heat transfer coefficient based on the out- 
side root tube area, A is the total heat transfer surface 
area of the exchanger used to compute U^, and AT^ is the 
mean temperature difference of the given exchanger. 

The object of the analysis, therefore, is to determine 

• • • 

Q 3 , Q 4 and given an initial listing of values for the 
design parameters. The list includes the following: 
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Tubeside mass flow rate, lbm/hr 

Entrance temperature of hot stream, °F 

Exit temperature of hot stream, °F 

Specific heat of hot fluid, BTU/lbm - °F 

Air mass flow rate, lbm/hr 

Entrance air temperature, °F 

Exit air temperature, °F 

Entrance air pressure, psi 

Specific heat of air, BTU/lbm-°F 

Tube inside diameter (ID) , in. 

Tube outside diameter (OD) , in. 

Fin height, in. 

Fin thickness, in. 

Fin spacing, in. 

Transverse Pitch, in. 

Longitudinal Pitch, in. 

Bank height, in. 

Bank width, in. 

Cross-flow arrangement 
Fin type 
Number of rows 
Number of passes 

Given heat transfer rate, BTU/hr 

Thermal Conductivity of tube material, BTU/hr .-ft .-°F 
Thermal Conductivity of fin material, BTU/hr .-ft .-°F 
Among the initial listing of design parameters above, 
there are parameters that are known and will remain constant 
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throughout the design problem. Also, there will be those 
parameters that are unknown and can vary, i.e. design varia- 
bles. On the way to determining the various heat transfer 
rates, other information will have been computed. This inform- 
ation includes the objective function, constraining functions, 
and other design data, such as the number of tubes per verti- 
cal row. 

The optimizer will then manipulate the design variables 
in order to find an optimum, while at the same time, performing 
a heat balance, that is: 

• * * • 

Q = Q 3 = Q 4 = Q 5 (15) 

where Q may be some given heat transfer rate. 

C. PERFORMANCE CALCULATION PROCEDURE 

With the temperatures, mass flow rates, and specific 
heats all specified in the listing of design parameters; 
whether they be constant or variable, the only unknown quan- 
tities on the right hand side of equations (12) through (14), 

are U , A, and AT . They will be determined as shown in 
m m 2 

Figure 1. 

1 . Mean Temperature Difference (MTD) 

For many flow arrangements, various approahces for 
determining MTD, mainly using diagrams, are available [36] , 
which have proven very useful in manual design efforts. For 
computerized design, however, an explicit, approximate equation 
is desirable in order to achieve a fast, sufficiently accurate 
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calculation of the mean temperature difference of a given 
flow arrangement. 

Roetzel, et al . [37], presented such an approximate 

equation together with empirical coefficients for nine 
counter-current cross-flow arrangements as they apply to 
air-cooled heat exchangers, 

Roetzel used the familiar equation for the MTD of the 
given flow arrangement fiT m : 



AT 



m 




(16) 



where AT, is the limiting case of pure countercurrent flow: 
m 



AT 



l 

m 





In 





) 

2 

) 



1 



(17) 



and F is a correction factor determined by a different set 
of coefficients for each flow arrangement. Roetzel reported 
the following function suitable for F: 



m n , 

F = 1 - £ £ a. . (l-v 1 ) sin(2 i arctan R) (18) 

i=l k=l X ' K m 



where is the dimensionless LMTD. 

m 



v 



m 




(19) 
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R 



( 20 ) 




T - T 
c 2 c x 

and the coefficients a. , of the approximating equation (18), 

1 , K 

were calculated using a standard least squares estimation 
program [38] , and are reported in reference [37] . The assump- 
tion that both streams were unmixed was used in their calcu- 
lation . 

When more than four tubeside passes are used, it is 
assumed that the heat exchanger has approached the limiting 
case of pure counterflow and F is set equal to one [2] . 

Having determined the MTD, the remainder of the analy- 
sis procedure follows Roetzel's [39] general approximation 
method for determining the mean overall heat transfer coeffi- 
cient, U^, for any flow arrangement while taking into account 
the variation of the heat transfer coefficients and the 
pressure drop with temperature and/or length of path. 

Before continuing with specific analysis procedures, 
a brief summary of Roetzel 1 s general approximation method is 
in order . 

The local overall heat transfer coefficient based on 
the outside root tube area can be written as follows: 



U = 



A " A In (r / r . ) T~ 

o 1 o o' 1 1 

A. H. 2irkL H n* 

11 o f 



( 21 ) 



where is the efficiency of the extended surface. 

In order to determine the individual convection heat 
transfer (film) coefficients, and H q , according to the 
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conventional methods, the coefficients would be considered 
constant, and the necessary fluid properties for their calcu- 
lation would be evaluated at some mean bulk temperatures, 

T, and T 
^b c b 

However, the film coefficients are not constant, but 

vary with temperature and/or length of flow path. Roetzel 

has taken these variations into account with the use of 

corrected reference temperatures. Two sets of corrected 

reference temperatures are determined: T, , T and 

n i C I 

T. , T . Therefore, for each set of corrected reference 

h n c n 

temperatures, the film coefficients are determined in the 
conventional manner using the reference temperatures in place 
of the bulk temperatures. 

With the film coefficients, H. , H , H. and H , 

X I °I 1 II °II 

two local overall heat transfer coefficients can be calcu- 
lated from equation (21), and U . 

Finally, the true mean overall heat transfer coeffi- 
cient is calculated as: 



tr - i [ ~ + — : i 

U m 2 Ui Un 



( 22 ) 



2 . Determination of Reference Temperatures 

The reference temperatures for a pure counterflow 
heat exchanger must first be determined from [39] : 



T' = T. + (T. -T. ) 

hj h 2 n^ h 2 



AT t - (T, -T ) 

1 h 2 C 1 
(T -T ) - (T. -T ) 
h l C 2 h 2 C 1 



( 23 ) 
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(24) 



T A 



II 



= T, + (T - T, ) 
h 2 h l h 2 



AT - (T -T ) 

II ^2 c ]_ 

( T -T ) - (T ~T ) 
h l c 2 h 2 1 



T ' 

c. 



T + (T - T ) 
C 1 C 2 1 



AT - (T -T ) 

1 h 2 C 1 
(T, -T ) - (T -T ) 
h l C 2 h 2 C 1 



(25) 



T' 

c 



II 



T + (T - T ) 
C 1 C 2 C 1 



AT tt - (T, - T ) 

11 h 2 C 1 

(T, -T ) - (T, -T ) 
h l C 2 h 2 1 



(26) 



where 



at = (T, -T )* 78868 • (T.-T )* 21132 

1 h i c 2 h 2 C 1 



(26a) 



AT 



II 



= (T, -T ) 



21132 



(T u -T ) 



.78868 



(26b) 



Equations (26a) and (26b) were derived by Roetzel in references 
[40] and [41] . For the special case where the fluid heat 
capacity rates, C, are equal, that is. 



m 



h p. 



= m 



the term 



AT. - (T, -T ) 

l h 2 C 1 
‘ Th 2 ' Tc 2 ) - <T h 2 - T = 1 > 



of equations (23) through (26) , with j 
.78868 or .21132 respectively. 



I and II, becomes 
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3 . Correction of Reference Temperatures for Given 
Arrangement 

With inlet and outlet temperatures fixed, pure counter- 
flow yields the highest mean temperature difference. There- 
fore, for any other arrangement, the temperature difference 
would be smaller. Thus, the corrections are applied in the 
following manner [39] : 



I 




I 



II 



T h 



- ip 



II 



h 



II 



T 

c 



I 




+ *c 



I 



'II 



T’ 

c 



+ <l» 



II 



II 



where \p, and ip 
y h, Y c 



are the temperature corrections of the 

k '"k 

hot and cold streams, respectively. The corrections are 



calculated as follows 



h j 



AT . 
3 



1 - AT /AT 
m x 



m 



1 + ‘W 



2/3 



c . 
3 



= AT 



1 - AT /AT 

m z 



m 



1 + (C c /C h ) 



2/3 



The corrected reference temperatures are now used to 
determine the thermal conductivities and absolute viscosities of 



38 



fluids for later use in the calculation of the film coeffi- 
cients. Thermal conductivity and viscosity data are usually 
presented in tabular form. However, for use on the computer, 
an explicit, simple approximate equation with temperature as 
the independent variable, was desirable. Water and air were 
chosen as two fluids that were likely to be involved in an 
air-cooled heat exchanger design. Figures 11, 12 and 13 
indicate that the thermal conductivities of air and water 
and the viscosity of air can be approximated by a second order 
polynomial. The viscosity of the hot tubeside fluid, water, 
must be treated specially, due to the following considerations. 

Calculations for the tubeside heat transfer coeffi- 
cient in the laminar, transition and turbulent regions all 

14 

involve the correction, (y/p ’ [42] . In the past, an 

iterative procedure was required to determine the inside tube 
wall temperature in order to evaluate • Roetzel [43], 

has developed an iteration-free method for determining this 
correction. Roetzel' s method assumes that the tubeside 
fluid's viscosity follows Andrade's Correlation, that is: 



Using the viscosity data of saturated water from reference 
[44] , the coefficients were determined through regression 
analysis with the resulting equation: 

M w = .01339 exp ( ^ ) (27) 
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where temperature is in degrees Rankine and viscosity is in 
lbm/ft-hr. Equation (27) was plotted next to the experimental 
data with results that indicate that water does indeed satisfy 
Andrade's equation, see Figure 14. 

4 . Uncorrected Tubeside Heat Transfer Coefficient 
The tubeside heat transfer coefficient is calcu- 
lated from one of three Nusselt-type empirical equations 
[42], as follows: 

For laminar flows, Reynolds Number < 2,100 



HD. 

l 

k 



1.86 Re 1 / 3 Pr 1 / 3 (^l 1 / 3 <— H ) ‘ 14 

L u wall 



(28) 



For transition regions, 2,100 < Re < 10,000 



HD . 

l 

k 



[Re 



2/3 



- 125] [1 + 



D. 

<x> 



2/3 



] Pr 



1/3 



.14 



wall 



(29) 



For turbulent flow, Re > 10,000 



HD . 
l 

k 



.023 Re* 8 Pr 1//3 




.14 



(30) 



Therefore, before any calculations can even begin, 
the Reynolds Number, Re, must be computed to determine the 
type of flow. 

The Reynolds Number will be calculated as: 




D i % 

A x u h. 

3 
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D i' m h is su PP lied in the initial list of parameters and 
is calculated from equation (27) . This leaves only the 
cross-sectional flow area, A , to be determined, where: 

X 



The number of tubes per pass, N , is determined geometrically, 

P 

having been given the tube bank height, number of rows, num- 
ber of passes, transverse pitch, diameters, and fin height 

initially. Figure 15 details the procedure of finding N . 

3 ? 

Mathematically , 



N . 1^22] th - D f - ^]/P t 

^ pass 



and would be a rational number. This rational number is 
then truncated to an integer number of tubes . 

With Re calculated, the Prandtl Number, 



Pr 



C P U h. 

k h. 

j 



is computed from given and previously determined thermo-physical 
properties . 

The uncorrected tubeside heat transfer coefficients, 

, that is; , without the factor (u/p w ) * ^ , can be 
calculated from the proper choice of equations (28) through 
(29) . 
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\ 




5. Calculation of Wall and Associated Resistances 



as : 



Equation (21) can be written in a more general form 



U 



R. + R , , + R 
1 wall o 



(31) 



where R. , R and R are the inside wall and outside heat 
1 w o 

transfer resistances respectively. Additional resistances, 
such as contact or fouling, can also be added here. 

In comparing equations (21) and (31) , the resistances 
can be computed as : 



R. 

l 



j 



^wall 




r 

o 



j 

In (D /D. ) 
o 1 



k wall 



Ola) 



(31b) 



R 

o 



j 




(31c) 



6 . Calculation of Airside Heat Transfer Coefficients 
In order to have the optimization program play a 
significant role in the selection of an optimized surface 
for a finned tube heat exchanger (which is one of the objec- 
tives of this project) , an explicit equation for H q involving 
tube and bank design parameters as independent variables is 
a necessity. 

In the past, comparison methods, as described by Shah 
[45] , were used to choose the best surface from a list of 
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surfaces for which experimental heat transfer and friction 
data existed. The data is presented in graphical form, 
where Friction and Colburn Factors are plotted as a function 
of Reynolds Number . With J, m^, c^ and Pr known, the film 
coefficient can be computed. 

Therefore, in previous heat exchanger optimization 
programs, a given surface described by its pitch, outside 
fin diameter, fin thickness, fin spacing, and outside tube 
diameter, had to be chosen beforehand. After the surface 
configuration had been selected, expressions for f and J 
were obtained by fitting polynomials to the experimental 
data described earlier. 

Briggs and Young [47] provide a means of getting 
past this obstacle with an improved convection heat transfer 
correlation for air flowing across triangular pitch banks of 
high finned tubes. 

Briggs and Young expanded the work of Ward and Young 
[48], to cover a total of 18 differently configured finned 
tube banks in order to determine the effect of fin thickness 
and tube pitch on the airside heat transfer coefficient. The 
heat transfer data for the high-fin tube banks were correlated 
to give : 

Nu = — 2 = .1378 Re* 718 Pr 1/3 (f-)* 296 (32) 

j 

where s is the distance between adjacent fins and £ is the fin 
height. Equation (32) is based on tubes having a wide range 
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of fin heights, fin thicknesses, fin spacing and outside 

tube diameter and can be used to predict H for a bank of 

o . 

3 

tubes six rows deep. Figure 11 of reference [48] is used to 
correct H' for banks of other than six rows. 

°j 

7 . Calculation of Fin and Surface Efficiencies 

In order to calculate the outside heat transfer 
resistance, which will be used to calculate the correction 
for the tubeside heat transfer coefficient and finally, the 
local overall heat transfer coefficients, the extended sur- 
face efficiency, n^/ must be computed. 

The surface efficiency accounts for the temperature 
drop from the root to the tip of the fin, due to the thermal 
resistance of the fin material. Thus, even though the heat 
transfer has been increased by the additional area of the 
extended surface, the area is not as effective as if it were 
to be all at the root temperature. 

The surface efficiency can be expressed as [49] : 

A f 

"* ■ 1 - x (1 -*' 

where <P , the fin efficiency for a radial fin, can be found 
from [ 49] : 



2 r 



m (r , 



- r 



I. (mr )k. (mr ) - k n (mr,)I.. (mr ) 

r 1 f 1 O 1 f 1 O i 

l I (mr )k. (mr,) + I. (mr,-)k (mr ) 
oolf ifoo 



The finned area, A^, and the total heat transfer area. A, are 
computed as follows: 
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A s = N T w n f i(D f 2 - D o 2 ) 



A = N T w[^-(D f 2 -D Q 2 ) + (l-n f t)D o „] 



With R calculated from equation (31c) , the correction to 
o 

the tubeside film coefficient can be made. 

8 . Correction of the Tubeside Heat Transfer Coefficient 
With the tubeside heat transfer coefficient dependent 
on the wall temperature, the dependence has either been 
neglected, or the wall temperature has been calculated with 
an iterative technique in the past. Roetzel [43] , has pro- 
posed an improved iterative-free method for finding the "Seider- 

14 

Tate" correction, (y/y , when the tubeside fluid 

wan 

obeys Andrade's viscosity correlation. 

From Roetzel 's formulation: 



where 



_y_ -) * 1 4 



wall j 



u . u . 

- -i + <-* 



+ V 



1/2 



u . 
3 



wall . 



,007 In ( 



v . = - . 07 In ( 




c . 

_L 



B .T. 

3 h. 



- 1 



c . 



V hj 



45 



where 



B . 
3 



h: 

i 



A (R 



wall 



+ R ) 



(Note that all temperatures are in degrees absolute.) 

All the parameters on the right-hand side of equation 
(21) are now available. Therefore, the two reference over- 
all heat transfer coefficients, and U^, can be calculated. 
The mean overall heat transfer coefficient follows easily 
from equation (22) . 

With U , Qc can be calculated, with the heat balance 
m ^ 

to be performed by the optimizer. 

9 . Calculation of Air and Tubeside Pressure Drops 

The final calculations before computing objective and 
constraint functions for the optimization problem involve 
the pressure drops in the heat exchanger. 

The basic equations that will be used for the calcu- 
lation of the pressure losses are as follows: 
a. Tubeside [42] 



A Pi 

where z 

2100. 



f . 

l 




w 





_u 

y 



z 



+ 



2 






(33) 



= .14 below Re = 2100 and z = .25 for Re greater than 



b. Air side [50] 




where A^ is the minimum flow area. 



( 34 ) 
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The friction factor for the tubeside flow, f^, is 
taken from Figure 9.5 of reference [42], assuming fully 
developed flow. For use on the computer, an explicit expression 
for f ^ was obtained by fitting a line and an exponential to 
the experimental data of Figure 9.5. This follows: 



Re <_ 1000 f. = .5/Re 

- 2694 

Re > 1000 f. = .003215 Re 

x 

Just as in the case of the airside film coefficient, 

for surface optimization on the computer, an explicit equation 

for the airside friction factor, f , is a necessity. Robinson 

o - 1 

and Briggs [50] , presented such an expression for f for air 
flowing across triangular pitch banks of finned tubes. Robinson 
and Briggs’ work closely parallels that of Briggs and Young 
[47]. The Robinson-Briggs Correlation: 



f 



o 



18.93 Re"* 316 (^-)~* 927 (—•) * 515 

o L 



(34a) 



covers the range of tube sizes and pitches used in air-cooled 
heat exchangers [50] . 

Therefore, with Re^ and Re.^, the four reference 
pressure drops, Ap. , and Ap , may be computed from equations 

lj °j 

(33) and (34). Following Roetzel's general approximation 
method : 
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For a gas, an additional correction is needed because the 
density in equation (33) is strongly dependent on pressure, 
which is changing through the exchanger. Using the inlet 
pressure as reference: 
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2Ap' , /9 

A p o = P 1 [1 - (1 - 1 (35) 



All the necessary information from an analysis view- 
point has now been calculated. Functions needed for the 
numerical optimization process shall follow. 

10 . Objective and Constraint Functions 

The objective functions available for minimization 
are defined as follows: 

a. Volume = wh[D £ + (n -1)P T cos 0] 

f r L 

where 



0 = arcsin (P /2P ) 

L. Li 



(36) 



(D 2 - D 2 ) tt . n , 

b. Heat Transfer Area = N T w[ j ■ ° + (1-— 'o 71 



Ap • m 

c. Air Horsepower = — — 

P 
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d. Airside Pressure Drop 

e. Tubeside Pressure Drop 

Limitations were imposed on the following inequality 
constraints in order to keep the design within practical 
physical bounds: 

a. The diameter ratio. 



DRATIO = jp- 
o 

must be kept reasonable. This can also be accomplished 
to some extent by placing side constraints on the 
design variables, £ and D q . See Figure 16. 

b. The optimizer must be prevented from driving the 
tube thickne s s , 



TUBTH = (D -D.)/2 
o 1 



to zero . 

c. The tubes must be kept from touching in both the 
longitudinal and transverse directions. 



TOUCHN = D f - P t 
TOUCHL = D f - P L 

TOUCHL and TOUCHN, therefore, must be kept below zero. 

d. Reasonable temperature profiles must be maintained 
at both ends of the exchanger. See Figure 10-7 of 
Ref . [44 ] , 



PROFH = 



T - T, 
c 2 h l 



PROFC = 



T - T. 
C 1 h 2 



that is, PROFH and PROFC, must be negative. 
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e. The number of tubes per vertical row, VROWR, shall 
not be allowed to go below 2. 

f. The free face area, that is, the minimum flow area 
for air, must obviously be greater than zero. 



DELSFF = [projected tube area] - hw 



where the projected tube area, STOTAL, is 

D f t w 

STOTAL = N [D w + — ] 

r o s 

and N is the number of tubes per vertical row. DELSFF 
must £e less than zero. 

g. The airside and tubeside pressure drops must be kept 
within design constraints. 

the argument of the square root of 



1 _ 2A P' 
p l 

has been driven below zero. It must therefore be 
constrained . 

i. To maintain an isoscele pitch bank, the angle, ©, 

as shown in Figure 17 and as defined in eq. (36) , must 
be constrained. One such constraining value is: 



h. From experience, 
equation (35) , 



ARG5 = 



© m = arcos (D^/2P l ) (37) 

The arguments of the arcs in and arcos of eq. (36) and 
(37) must be constrained from going beyond 1. 

j. The heat balance described by eq. (15) is performed 
by three equality constraints: 



• • 

Q RATIO = Q 5 /Q (38) 

QRATI01 = Q 3 /Q 5 (39) 
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QRATI02 



(40) 



= q 4 /q 3 

All constraints are set equal to 1. 

Besides the constraints described in paragraphs 
j . , above, lower and upper bounds are placed on the 
variables to assure a reasonable design. 



a . through 
design 
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IV. RESULTS 



A. BACKGROUND 



Case studies were chosen as the best way to test the 
capabilities of the program for Heat Exchanger Design using 
Numerical Optimization (HEDSUP) . The design problems posed 
were made as realistic as possible. 

1 . Capabilities 

HEDSUP currently has the capability to design for 
nine different configurations of triangular pitch banks of 
finned tubes: 



TYPE 


1 


— 


1 


ROW, 


1 


PASS 


TYPE 


2 


- 


2 


ROW, 


1 


PASS 


TYPE 


3 


- 


3 


ROW, 


1 


PASS 


TYPE 


4 


- 


4 


ROW, 


1 


PASS 


TYPE 


5 


- 


2 


ROW, 


2 


PASS 


TYPE 


6 


- 


3 


ROW, 


3 


PASS 


TYPE 


7 


- 


4 


ROW, 


2 


PASS 


TYPE 


8 


- 


4 


ROW, 


2 


PASS 


TYPE 


10 


— 


PURE COUNTERF 



TYPE 10 will include exchangers with a configuration of n rows, 
n passes, where n can go from five to 20. 

The banks must be constructed of high-finned tubes 
U > .0625 in. [42]) with the fins having a rectangular pro- 
file of constant thickness. Additional profiles can be inserted 
into HEDSUP quite simply, provided that its fin efficiency can 



52 



be expressed explicitly as: 



n = f(£,t,H,k) 

See Subroutine FINEFF of the program listing, Appendix D. 

At present, HEDSUP can provide the design parameters 
for an air-cooled heat exchanger optimized for any one of the 
following design objectives: 

(a.) Minimum Volume 

(b.) Minimum Heat Transfer Surface Area 
(c.) Minimum Air Horsepower 
(d.) Minimum Airside Pressure Drop 
(e.) Minimum Tubeside Pressure Drop 
Additional design objectives can be used, provided that they 
can be expressed explicitly as a function of the design varia- 
bles and they are added to the common block. It should also 
be pointed out that any design variable may simultaneously be 
a design objective as long as it conforms to the restrictions 
of both. For example, an exchanger may be designed for mini- 
mum bank height. 

The airside fluid is restricted to dry air. The tube- 
side fluid is presently limited to water in single phase. 

Other tubeside fluids can be used by HEDSUP, provided that 
their visocisites obey Andrade's Law and the fluid's thermal 
conductivities can be expressed explicitly as a function of 
temperature. The fluid's specific gravity would also have to 
be placed in the denominator of eg. (33) . 
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B. CASE STUDIES 



1 . Case One 

a. Problem Formulation 

An air-cooled heat exchanger is to be designed 
for minimum volume with a heat transfer rate of 10,000,000 
BTU's per hour. Water is to be coaled from 200 °F to 125 °F. 
Dry air will enter the exchanger at 95°F and leave at 130°F. 
Specifications call for a fan that can produce a pressure 
difference of two inches of water. 

b. Design Variable Framework 

From the list of design parameters in Section 
III.B, the design variables must be singled out, identified 
to COPES, and given side constraints. All parameters must 
be given an initial value. Only the values of the design 
variables will change. 

Assuming constant specific heats. 



c =1.0 BTU/lbm- °F 
Pw 




.24 BTU/lbm- °F 



the mass flow rates of both fluids can be determined from 
eqs. (12) and (13) , because the heat transfer rate and tern' 
perature differences are given. 



m 



w 



m 



a 




133333. lbm/hr 



190476.2 lbm/hr 
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A cross-flow arrangement, fin profile, tube 
material, and fin material must be chosen. 

The design variables for this example are, there- 



fore : 



.232 < D. < 2.325 in. D. 1 

l l 

.24 < D < 2.5 in. D 1 

o o 

.0625 in. < i < 00 2. 1 

.01 < t < .0235 in. t 1 

.08 <S <.125 in. S 1 

0.0 < P <4.0 in. P 1 

Li L 

0.0 < P fc < 4 .0 in. P t 1 

0.0 < w < 500 in. w 1 

0.0 < h < 500 in. h 1 



2.0 inches 
2.5 inches 
.46 inches 
.023 inches 
.111 inches 
2.125 inches 

4 . 0 inches 
490 inches 
350 inches 



The side constraints on the design variables are of a practi- 
cal nature with the exception of the lower bounds on fin 
height. Recall that the use of eq. (32) is restricted to 
high-fins. High-fins will also tend to keep the fluid unmixed, 
which was an assumption used when defining the coefficients 
a. ,, used in eq. (18) . 

J. , K 

c. Constraint Framework 

From the problem statement, the air side pressure 
drop must be less than two inches of water or .0722 psi. 
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0 < 


Ap < 

*a 


0722 psi 


From a 


practical 


standpoint: 


O 

• 

o 


< 0 < 1 


.3 


1.0 


< DRATIO 


< 2.5 


.018 


< TUBTH 


< .18 inches 


— 00 


< p < 

^w 


.14 psi 


The equality constraint 




°5 

Q 


1.0 


°5 = 


= U A AT 
m 





Q = 10,000,000 Btu/hr 

will satisfy the heat balance. 

The bounds on TOUCHN, TOUCHL ARG5 , ARG7 , ARG8 , 
DELSFF and VROWR were discussed in Section III. C. 10. Con- 
straints on PROFH and PROFC were unnecessary because all 
temperatures were specified in the problem statement. 

As with design variables, constraints must be 
identified to COPES by location in the common block. See 
the global catalog. Appendix A. 
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d. Methodology 

Ideally, a three dimensional design matrix can 
now be constructed of optimum exchanger designs with minimum 
volumes. The matrix would be constructed by first holding 
the tube and fin materials constant and varying the configura- 
tion, i.e. Type 1, Type 2, Type 3, etc. Next the tube material 
would be varied with the fin material and exchanger configura- 
tion held constant and so forth. However, for case study one, 
the tube material will be chosen as copper, k = 200 BTU/ft-hr- °F , 
and the fin material will be aluminum, k = 118 BTU/ft-hr- °F . 

Also, in order to simulate an actual trade off 
study, the constraint framework will be fixed throughout the 
individual case studies. 

Problems arise in constructing the matrix when 
trying to determine the true minimum volume design for each 
configuration. Unfortunately, the choice of initial design 
parameters (starting point) , coupled with the input parameters 
for ALMM, will sometimes lead to entirely different optimum 
designs with volumes differing by over 100%. 

The ALMM parameters include the initial multiplier, 
CC, the multiplication factor, CMULT, and the maximum multi- 
plier value, CCMAX. Experience has shown that setting 

CMULT = 2.0 

CCMAX = 1000. 

will suffice for almost all problems. However, there is much 
"artwork" involved with the choice of CC . From experience. 
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an initial multiplier of 10 works well when starting far 
from the equality constraint, i.e.: 

.8 < Q RATIO < 1.2 

However, when approaching the heat balance, i.e., the equality 
constraint, a CC = 100 is necessary for the heat balance to 
converge. Equation (15) is considered satisfied when the heat 
transfer ratios (eqs. (38) through (40)) are less than 0.1 
percent. 

Therefore, it is obvious why the chosen initial 
design is so critical. Together with the choice of CC it 
will determine how and to what design the optimizer will con- 
verge. As an example, see Table 1. Notice the calculated 
heat transfer for the initial design. This value is the product 
of U m , A, and AT^ calculated by ANALIZ using the initial design 
parameters, some of which are mere estimates. Recall it is 

the job of the optimizer to vary the design variables in order 

• • • • 

to bring Q<- equal to Q, and Q 4 and at the same time minimize 
the objective function. 

In order to remove some of the "artwork" and 
try to ensure a true optimum design, the following procedure 
is recommended to generate the design matrix: 

(1) Begin with a Type 2 configuration; input the initial 
design and constraint values enumerated in sections 
(a) through (d) above; let CC = 10.; execute. 

(2) If the heat balance of the resulting design has not 
converged, but is within 20%, use the design results 
as a new starting point and set CC = 100 (if the heat 
balance is not with 20%, let CC = 10.). 
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TYPE 7 EXCHANGER CONFIGURATION OPTIMUM DESIGN 





RUN 


I 


RUN II 


DESIGN 

VARIABLES 


INITIAL 

DESIGN 


OPTIMUM 

DESIGN 


INITIAL 

DESIGN 


OPTIMUM 

DESIGN 


D . , inches 

l 


. 6786 


.5392 


2.0 


1.168 


D q , inches 


.7201 


.5767 


2.5 


1.278 


SL , inches 


. 1618 


. 128 


.46 


.294 


t , inches 


.0217 


.0208 


. 023 


. 023 


S , inches 


.08 


. 08 


. Ill 


. 08 


P^ , inches 


1.044 


1.02 


4.0 


2.557 


P , inches 

Li 


1.044 


.841 


2.125 


1.899 


h , inches 


342.1 


233.4 


350. 


172.9 


w , inches 


286.6 


186.2 


490. 


253. 












Q , BTU/hr 


21,988,624. 


9,999,925. 


20,650,352. 


9,999,811. 


3 

Volume, ft 


213.12 


72.92 


553.2 


153.83 



Table 1 
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(3) Repeat step (2) until convergence. NOTE: CC may 

be adjusted up to 150 when approaching convergence. 

The design should converge following the use of 

CC = 100. If too much adjustment of CC is necessary, 
reaching the optimum from the starting point is 
unlikely . 

(4) Ensure a minimum design by beginning step (1) with a 
different initial design. 

(5) After finding an optimum design for a Type 2 con- 
figuration, use that design for the starting point 
for a Type 3 configuration. This assures a reason- 
able starting point, probably close to the optimum 
for Type 3 . 

(6) Repeat steps (1) through (3) for a Type 3 configuration. 

(7) Repeat steps (5) and (6) for the remaining configuration 

e. Design Matrix 

The design matrix is presented in Table 2. The 
optimum design is a Type Four configuration and is shown in 
Figure 18. Typically, for this case study, when starting far 
from the final design, COPES/CONMIN would require approxi- 
mately 1700 calls to ANALIZ to reach an optimum. However 
when beginning from a reasonable starting point with CC = 100, 
only 600 calls were needed. Note that each call to ANALIZ 
requires approximately .06 seconds of CPU time on an IBM 360/67 
2 . Case Study Two 



a. Problem Formulation 

An air heater is to be designed to fit into a 
space 8' X 24' X 4'. The heat exchanger is to heat 1,000,000 
lbm/hr of dry air from 75°F to 130°F. 256,000 lbm/hr of water 

at 200 °F is available. Design the heater so that the required 
air horsepower is at a minimum. 
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b. Design Variable Framework 



Assuming constant specific heats. 



c =1.0 BTU/lbm- °F 

Pw 



c = .24 BTU/lbm- °F 

P a 



eqs . (12) and (13) will yield the required heat transfer 

rate and the outlet water temperature. 



Q = mC (T -T) = 13.2x10 BTU/hr 

p a c 2 1 



“h C p T h. ' Q 

MV 1 



= 148.44 °F 



m. 



Therefore, the design variables are as follows: 



232 < D. < 2.325 in. 

l 



D i = 2.0 in. 



25 < D < 2.5 in. 

o 



D = 2 . 5 in . 

o 



.0625 


in . < 


£ 


< 00 




.46 in. 


i — 1 

o 

• 


< 


t 


< 


.0235 in. 


t 1 = 


.023 in. 


• 

o 

00 


< 


s 


< 


.125 in. 


s 1 = 


.111 in. 


o 

• 

o 


< 


P L 


< 


00 


V ■ 


2.125 in 


o 

• 

o 


< 


P t 


< 


00 


p t i ■ 


4.00 in . 
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c. Constraint Framework 



When designing for minimum horsepower, the opti- 
mizer will naturally try to drive the design to a maximum 
volume in order to reduce airside pressure losses. Therefore, 
it is reasonable to constrain volume as follows: 

0 < Volume < 768 ft^ 

The other constraints are: 

0.0 < 0 < 1.3 



1.0 < DRATIO < 2.5 



.018 < TUBTH < .18 in. 

0 < Ap^ < .14 psi 



and for the heat balance: 




1 



d. Design Matrix 

With the tube and fin materials fixed, as in case 
one, the matrix is presented in Table 3. The optimum design 
is a Type 4 configuration, as shown in Figure 19. In case 
study two when starting far from the optimum, the optimizer 
called ANALIZ approximately 1900 times. When starting close 
to the final design, for example using the design for a Type 3 
configuration as a starting point for the Type 4 design, 
COPES/CONMIN only required 596 calls. 
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Table 3. CASE II DESIGN MATRIX 
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V. CONCLUSIONS 



The intent of this investigation was to couple an analysis 
program with a numerical optimization scheme, COPES/CONMIN, 
to produce a complete, detailed design program for an air- 
cooled heat exchanger, (HEDSUP) . In addition, the analysis 
program was to be written such that: 

1. the variation in the film coefficients with temperature/ 
length of flow path would be taken into account 

2. the surface would be optimized 

3 . it would be iterative free and thus minimize the CPU 
time required during an actual trade off study. 

The results from test cases using ANALIZ coupled with 
COPES/CONMIN in its present form were unsatisfactory. Although 
COPES/CONMIN could optimize the objective function satisfying 
the inequality constraints, a reliable heat balance could not 
be obtained. The solution to this problem was the addition 
of the ALMM option to COPES. In this way, the method of 
feasible directions, which works best with inequality con- 
straints, was used to satisfy the inequality constraints. 

The multiplier method, which works best with equality con- 
straints, was used to perform the heat balance. 

The results of the case studies show that HEDSUP will 
yield reliable designs for various design objectives and 
problems with only some trial and error application of the 
initial Lagrange multiplier. However, precautions must be 
taken to overcome the relative minima that plague this design 
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problem. Table 1 shows vividly the problem of relative 
minima with one "optimum" design having a volume over 100% 
greater than the "true" optimum. 

The value of numerical optimization in a design problem 
of this size cannot be overemphasized. For example, in case 
study 1, the problem is taking place in a nine-dimensional 
design space and intuition on how an optimized design "should" 
turn out is quickly lost. Figure 20 helps to illustrate this 
point. When only varying two design variables, h and w, 

Figure 20 shows that, when beginning from the initial starting 
point with a Type 7 configuration, a design satisfying both 
the equality constraint and the air pressure drop constraint 
could never be found. However, Table 2 indicates how the 
numerical optimization routine has varied the other seven 
design variables in order to shift the constraints and yield 
an optimized design. 

The results show that AHDOP did vary the surface design 
variables: D^, D q , l , t, S, P fc and P^, in order to produce 
an optimum heat exchanger. This capability of surface opti- 
mization is dependent upon the use of the Briggs-Young and 
Robinson-Briggs Correlations, eqs . (32) and (34a), respec- 

tively. The reliability of the correlations as compared to 
the "conventional" method is questionable. Actual experimental 
data for a particular tube and pitch will always be the most 
useful in predicting pressure drop and film coefficients of 
across banks of finned tubes. However, the correlations 
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mentioned above cover the ranges and pitches used in air- 
cooled heat exchangers, and should therefore be sufficiently 
accurate in predicting H q and Ap a « 
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VI . RECOMMENDATIONS 



In addition to the insight that this investigation has 
given into the generation of an automated air-cooled heat 
exchanger design, it has also generated an awareness of this 
investigation's shortcomings. As mentioned in the review of 
previous work in this area, each of the optimization methods 
has its own limitations; none is completely general. Pre- 
sented herein are recommendations for improving upon and 
furthering development of HEDSUP . 

1. HEDSUP should be expanded to include the capability 

for two-phase tubeside fluids. Mott, et al . [2] discusses a 

method involving two-phase tubeside fluids that would be 
compatible with HEDSUP. The modular design of ANALIZ will 
aid this effort. 

2. Research in the area of numerical optimization using 
discrete variables would benefit HEDSUP immensely. With dis- 
crete variables, the design of the exchanger could be accom- 
plished with "off-the-shelf" materials. At present, the use 
of the optimizer is restricted to continuous variables with 
continuous first derivatives. The ability to work with dis- 
crete variables would also eliminate the need for the design 
matrix. The optimizer could optimize for type configuration, 
fin profile, fin material and tube material. 

3 . Additional research with ALMM is needed to remove the 
"artwork" involved with choosing an initial multiplier and 



68 



thus increase reliability and hence reduce CPU time. The 
research should be concentrated in two areas: 1) debugging 

the new optimizer with AFMM and 2) scaling of design variables. 

4. The addition of cost as an objective function would 
increase the attractiveness of HEDSUP. Mott, et al. [2], 
and Fontein and Wassink, [18] , have presented much useful 
information in this regard. 

5. Mechanical constraints such as tube bursting stress 
and tube vibrations should be included in HEDSUP. 
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VII. FIGURES 




Figure 1 
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Figure 10. Configuration of Air-Cooled Heat 
Exchanger 
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BflMK WIDTH 



APPENDIX A 



HEAT EXCHANGER DESIGN USING NUMERICAL OPTIMIZATION (HEDSUP) 

USERS MANUAL 

In order to execute HEDSUP, it is necessary to provide 
formatted data for COPES [25] , followed by formatted data 
for the ANALIZ portion of the program. This section defines 
the data which must be supplied. The data is submitted in 
segmented blocks. All formats are alphanumeric for TITLE 
and END cards, F10 for real, and 110 for integer data when 
supplying COPES. For ANALIZ, formats are F14 for real and 110 
for integer data. 

While the COPES user's sheets define the data in formatted 
fields of ten, the COPES manual does provide means of simpli- 
fying this task through unformatted data input. 

The included GLOBAL CATALOG defines objective functions, 
design variables and constraints along with their location 
for ease in compiling the necessary COPES data. 
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GLOBAL CATALOG 



GLO BAL 
LOCATION 


FORTRAN 

NAME 


DEFINITION 


1 


FMDOT 


Fluid Mass Flow Rate, lbm/hr 


2 


FLTMP1 


Inlet Fluid Temperature, °F 


3 


FLTMP2 


Outlet Fluid Temperature, °F 


4 


AMDOT 


Air Mass Flow Rate, lbm/hr 


5 


ARTMP1 


Inlet Air Temperature, °F 


6 


ARTMP2 


Outlet Air Temperature, °F 


' 7 


TUBEID 


Tube Inside Diamter, in 


8 


TUBEOD 


Tube Outside Diameter, in 


9 


FINHT 


Fin Height, in 


10 


FINTH 


Fin Thickness, in 


11 


FINSP 


Fin Spacing, in 


12 


PITCHL 


Longitudinal Pitch, in 


13 


PITCHN 


Transverse Pitch, in 


14 


BANKW 


Bank Width, in 


15 


BANKH 


Bank Height, in • 


16 


VOLUME 


Volume, ft~* 


17 


AREA 


o 

Heat Transfer Area, ft 


18 


FHP 


Air Horsepower, HP 


19 


PPA 


Airside Pressure Drop, psi 


20 


PPW 


Tubeside Pressure Drop, psi 


21 


DELQl 


Q - Q 5 , BTU/hr 


22 


DELQ2 


Q 5 - Q 4 , BTU/hr 


23 


DELQ3 


Q 4 - Q 3 , BTU/hr 


24 


DRATIO 


D f/ D o 


25 


TUBTH 


Tube Wall Thickness, in 


26 


TOUCHN 


Tip-to-Tip Clearance, Transverse 


27 


TOUCHL 


Tip-to-Tip Clearance, Longitudinal 
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GLOBAL CATALOG (Cont.) 



28 


PROFH 


Tc n - Tin > ° F 


29 


PROFC 


T cl -T ho , °F 


30 


ARG5 


1 - [2Apq/poo] 


31 


QDOT1 


Q s , BTU/hr 


32 


DELPHI 


0 -- e raax’ radians 


33 


QRATIO 


q 5 /q 


' 34 ‘ 


ARG7 


P r /2P L 


35 


ARG8 


D f /2P L 


36 


VROWR 


Number of Tubes per Vertical Row 


37 


DELSFF 


Negative of Free Face Area 


38 


THETA ■ j 


| arcsin (ARG7) 














































































* 
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DESCRIPTION: Title card 
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nnsCRIPTION: Program Control Parameters. 
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DESCRIPTION: Integer optimization control parameters. 
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constraint values during the one-dimensional search. 



UMAX: Maximum number of optimization iterations allowed. DEFAULT = 20. 

ICND1R: Conjugate direction restart parameter. DEFAULT = NDV + 1. 

NSCAL: Scaling parameter. GT.O - Scale design variables to order of magnitude 
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1) Currently NFDG must be zero in COPFS. 



PE.SCK1 1*T ION; Floating point optimization program parameters. 
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FDCUM: Minimum absolute step in finite difference gradient calculations. 

DEFAULT = 0.001. 



CT: Constraint thickness parameter, DEFAULT « -0.05. 

CTMIN: Minimum absolute value of CT considered in the optimization process. 

DEFAULT « 0.004. 
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1) The DEFAULT values for these parameters usually work well. 



DESCRIPTION : Total number of design variables, design objective identification and sign.- 
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+1.0 indicates maximization. -1.0 indicates minimization. If SGNOPT is 



DESCRIPTION: Design variable bounds, initial values and scaling factors. 
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DESCRIPTION : Design variable identification 
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DESCRIPTION : Number of constrained parameters 
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DESCRIPTION: Constraint: identification and constraint bounds. 
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linear constraints , LCON - 2 indicates equality constraint 



bound on the constrained variables. If BL. LT . -1. OE+15, no lower bound 
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DESCRIPTION : COPES data 'END* card 










CO 








«H 








J= 








4J 






• 


CO 






6? 


5 






4J 


o 






3 


rH 






•a 


*H 








O 






CO 


CM 






Cd 








cu 


A 


CO 




o 


SI 


H 




a 


M 


2 








U4 




<u 


< 


H 




<— 


2 


2 




u 


< 


C 


• 






o 


cn 




A 




i 


O 


3 




pH 




c 






*3 


•H 




CO 


C 


« JJ 




c 


3 


3 3 




a 




-u o 




a 


a 


S3 U 




iH 


xz 


*3 




o 


u 


"O 




a 




<U 0) 






u 


3 *H 




c 


3 


CU iH 




■H 




c CU 






u 


«H C. 




• 


<3 


3 




a 


<0 


CO 3 




2 


a 


1 




3 


Cm 


^4 W 






3 


o 3 








CJ CO 




•3 




3 




W 


CO 


3 




O 




XZ 3 








M ^ 








M 




a 


-a 


31 






j- 


*3 M 




H 


(3 


c o 






U 


3 <M 






0) 


CO 3 




CO 


*H 


*M M 




uc 


XI 


JT 3 


c 


erf 


H 


c- a 




< 






1U 


-H 2 








SsJ 


/*N 


y— N ^“N 




2 


H 


cn m 
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ANALIZ DATA 



Cfl 

5-1 

0) 

*J 

0) 

E 

3 

5-i 

03 

fa 

a 

•H 

Cfl 

a 

.o 

3 

H 



2 

O 



H 

PUi 





}-i 




fa 






fa 


0 






0 






B 




a 


fa 


rO 


•s 


0) 


0 


1 1 


0) 


5-i 


1 




5-1 


3 


E 


•S 


3 


4-> 


.a 


<U 


u 


3 


r— ! 


u 


03 


5-i 


\ 


03 


5-1 


a 


!3> 


5-1 


a 


fa 


H 




fa 


E 


fa 


3 


E 


a 




O 


OJ 


•u 


r\ 


r— 1 


*_) 




u 


4-1 




03 


3 




03 


*H 


a 


w 


*H 


3 


x 


w 


3 


r— i 




03 


r— I 


4-1 


o 


E 


4-1 










u 


4-1 


03 


U 


a) 


*H 


•H 


<U 


r— 1 


a 


3 


r— 1 


•U 


a 


i — 1 


c 


3 


a 


4-^ 




O 


w 


1 


r— I 


1 

CNJ 


i 


• e~ 


X 

H 


X 

H 


p_ 

CJ 



C/3 



fa 

fa 
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1) This begins ANALIZ data input. 

2) Input must be formatted. 

3) Comment cards are no longer permitted. 



DESCRIPTION: Airside Parameters 



o 

Pm 





5F14.4 










70 


3 

a 




LT, 


03 

CL 

CD 




CNJ 

<r 


X 




28 


rH 

a 

H 




<r 

i-H 


m 

a 





Li 




Pm 






DC 


Pm 


o 


Pm 






o 




o 




6 




* 


I 




42 


r 


qj 


2 




«— 1 


0) 


Li 




•H 




Li 


0 


rH 


CO 


•s 


0 


mi 




PL 


QJ 


M> 


cd 


3 




mi 


cd 


Li 


H 




c0 


Li 


QJ 


42 


QJ 


L 


QJ 


CL 




L 




CL 


2 


* 


D 




2 


QJ 


M 


CO 


o 


a 


Ml 


cd 


CO 


r—\ 


u 




QJ 


aj 


CM 




L 


42 


u 




mi 


•rH 




CL 


co 


QJ 


cd 


a 




co 


T-H 




•H 


Ml 


cd 


c 


Ml 


CM 


C 


2 


•H 


QJ 


•rH 


QJ 






r-H 


a 


•H 


Li 


Li 


U 


QJ 


42 


■H 




D 


CL 


2 


cd 


cd 


O 


co 


cd 



III!! 

.H CNI Cd 

co a a a. 8 

£2 fH E— i CJ cl 
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I 



DESCRIPTION: Tube Geometry 
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DESCRIPTION: Tube Arrangement 





4F14.4 










vO 

LO 


5 




42 






28 


hJ 

PH 




14 


4-1 

Ol 







CD 








CU 






C/D 


X 






CU 


a 






X 


a 






a 


♦H 






c 








•H 




CD 








CU 


CD 


0S 


CJ 


JC 


CU 


r* t 


4J 


o 


x 


a 


•rH 


d 


CJ 


U 


a- 


•H 


d 


•H 






•H 


Cu 


r— 1 


0S 






d 


4-4 


*\ 


cu 


d 


f 




CD 


♦H 


'So 


4-4 


l-i 


T3 


•H 


T3 


CU 


d 


cu 


•H 


> 


4-4 


pd 




CD 


•H 






d 


CO 






a 


d 


d 


d 




o 


03 


03 


4-4 


t— i 







Q 

hJ 

W 

M 



1111 
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DESCRIPTION : Integer Parameters 



o 

rH 

M 



LT| 



CN 



00 

CN 



co 

CO 

< 

pH. 

2 : 



CO 

3 = 

O 

04 

Z 



W 

P-. 

>* 

H 



W 

CL. 

H 



Q 

hJ 

W 



G 

O 

•H 

4J 

CD 

U 

a 

00 

•H 

M— 1 

G 

O 

o 



w 

pH. 

>• 





CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CD 


CD 


CD 


CD 


G- 


G, 


C- 


G* 


rH 


rH 


rH 


rH 






* 




CO 


CO 


CO 




£ 


s 


2 


0 


0 


0 


O 


u 


u 


jh 


u 


r—i 


CN 


cn 


<r 



0} 

a) 

cn 

co 

CD 

C. 



CO 

o 

JH 



CO 

2 

o 

u 



CO 

o 

5-4 



CN co 



I I I I I I I 

r— I CNi CO <f LO VO r^* 



109 



3 passes 



crossflow 



CO 

Q) 



o 



2 
o 
1—1 
< 4-4 

u 

Q) 
u 
C 
3 
O 
O 

<u aj 



u 

03 

pH 

3 

CO 

C 

03 



Q 

W 



CO 

a> 





u 


U 




a 


CO 


CO 




3 


3 




aj 




CO 


<r 


Cu 


a. 




u 


o 


03 








Q) 




u 


a 








rH 








I 


1 


i 


*H 


i 


<4-4 


<4-1 








<4H 




o 


o 






o 


o 








co 


ON 


i— i 


JmI 


i— i 












a, 




<1) 


aj 












r3 


& 








d 




B 


G 








iH 




3 


3 








<4-1 

I 




C 

i 


d 

1 








1 

w 




i 

CO 


1 

CO 








PU 




3 


CO 












o 


< 








H 




PS 


PU 








■""D 




2 


z 
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DESCRIPTION: Miscellaneous Parameters 



H 



O 

P^ 



rH 

pH 

<r 








3 




N 

d- 

2 
i— i 
Ph 




30 

N 

w 

XI 

O 

H 




<r 

H 

•O' 







pH 








0 


PH 






1 


o 






p 


1 






f-* 


p 






1 








u 


1 






CM 


J-J 






— 


CM 






O 








H 


o 






X 


H 








X 






#* 




ex 






#* 


o 




cd 


1 1 


I- 




•H 


cd 


£ 




P 




X 




01 


p 


i-H 




4-> 


0) 


— ■**. 




cd 


4-J 


CM 


P 


e 


2 


X> 


jo 




£ 


1 — 1 


\ 


a; 




1 


o 


x> 




4-4 


H 


D 


•H 


CM 


XI 


ij 


Mm 










<r 




CM 


CM 


ro 


OJ 


o 


o 


• 


u 






m 


cd 








Pi 


u 


4-4 






•H 


•H 




P 


> 


> 


p 


a; 


•H 


•rl 


•H 


cm 


4-J 


4-4 


cd 


c n 


o 


a 




£ 


3 


£ 


CM 


cd 




73 


o 


p 


c 


£ 




u 


o 


O 


4-4 




a 


O 


£ 


u 






cd 


cd 


i-H 


i— 1 


4-4 


<u 


cd 


cd 


cn 


x: 


e 


B 


£ 




p 


U 


O 


£ 


0) 


V 


a 


<D 


x: 


x: 




> 


4J 


4-J 


w 


•H 






cd 


to 


1 


1 






w 






1 


X 


z 


I 




o 


M 






H 


Ph 




•O' 






cX 
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DESCRIPTION; Penalty Parameters 
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DESCRIPTION : Optimization Parameters 






H 



O 

fa 




o 

II 

>-< 

li 

>-• 

li 

on 

>-< 

II 

CN 

JH 

II 
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If NCALC = 3 let Y1 



DESCRIPTION: Augmented Lagrangian Function Method Data 
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COPES DATA 



APPENDIX B 



SAMPLE USER'S, INPUT 
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DATA BLOCK D - OMIT IF NDV 





O 

II 

> 



0 

1 



u* 






CJ 

o 

S3 

< 



e- 




H 

s 

s 

*-* 

o 

CJ 


§ 

s 

o 

u. 


o 

<— i 
►*< 
vT 




'/> 








SCAL 


O 

o 


O 

o 




2.0 


1 2.5 


VUB 


2.325 


2.5 


VLB 


.232 


.25 
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DATA BLOCK. F - CONT. 
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DATA Ul.OCK G - OMIT IF NDV 



C-t 




o 


z: 


H 


fH 


w 


d 


U. 




S 


o 


o 


o 


.-H 


o 


u< 








CN 



o 

II 



ZD 

< 



O 

o 



o 

o 



o 

o 



o 

o 



o 

o 



o 

o 



o 

o 



o 

o 



o 

o 



o 

oo 

Q 



CO 



<7\ 



CM 



<T) 






r: 

o 

00 



CN 



cn 



m 



vO 



co 



cr\ 
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DATA CLOCK H - OMIT IF NDV 




119 



DATA BLOCK I - CONT 

































































SCAL2 


o 

d 










SCAL2 


O 

o 










SCAL2 


o 

o 




2 


00 




LC0N 


O 




2 


o 

o 




LC0N 


o 




3 

CQ 


.0722 




SCALl 


O 

o 




JC0N 


CN 




SCALl 


o 

• 

o 




o 

CJ 

*-> 


r~! 




SCALl 


O 

• 

o 


CO- 




00 

O 


co- 


ICOI'! 


kO 

CN 


</> 


CO 


-1.0 + 20 


< 0 - 


8 

'O 


ON 


<ry 




-1.0 + 20 
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DATA BLOCK I - CONT 




121 



DATA BLOCK I - CONT 




122 



DATA BLOCK I - CONT 




12 3 



DATA DI.OCK V - COPliS END OF DATA CARD 







% 
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ANALIZ DATA - CONT 



H 

<r- 

£ 

O 

lu 


3F14.4 


FORMAT 


5F14.4 


FORMAT 


3F10 . 0 










H 

3 

2 

o 

a- 










































' 








































O 

O 






















<r 

>• 


O 

O 


















R3 


O 

O 


on 

>* 


o 

o 


CCMAX 


1000. 














R2 


O 

O 


CM 


o 

o 


H 

hJ 


O 

CM 






















O 
















R1 


O 

o 


•H 


O’O 


a 

a 


10.0 
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APPENDIX C 



SAMPLE OUTPUT FROM COPES 



AAAA AAAAAAAA AAAAAAAAA AAAAAAAAAAA A A A A A A A A A A A A. A A AAA A AA A A A A A A A A AAA A A A AAAAAAAAAAA/ 



03/16/flO 13 . 5*« . > l 



FILE: SAV 
1 



DATA 



T 1 



NAVAL POSTGRADUATE SCHOOL 



ccccccc • 

C 

£ 

c 

c 

ccccccc 



OOC.OCOO 
0 0 

8 8 

L) 0 

o a 

OGCCCOQ 



ppppppo 

p p 

p p 

ppppppp 

p 

p 

p 



EEEEEEE 
E . 

£eee 

E 

E 

EEEEEEE 



sssssss 

s 

Issssss 

s 

s 

sssssss 



CONTROL PROGRAM 
FOR 

ENGINEERING SYNTHESIS 



TITLE 



AIR 


HEATER DESIGN - OPTIMIZE 


VOLUME . 


CARO 


IMAGES CF CONTROL DA FA 




CARO 


IMAGE 


. 


1 ) 


AIR HEATER DESIGN - OPTIMIZE VOLUME. 


2) 


2,9 




2 ) 


2 9 




3) 


* BLOCK C 




A) 


Of 8 C f Of 10 




4) 


0 30 


o 

o 


5 ) 


$ 3LGCK 0 




6) 


0. 




6) 


0. 




7) 


0. 




7) 


0. 




8 ) 


$ BLOCK E 




9) 


9 f l 6 f — 1 • 




9) 


9 16 


-1 . 


10) 


S BLOCK F 




11 l 


* 232 t 2 • 325 f.6736 




fit 


. ^ ^32 2.323 

•23f2.3f.7201 


.6786 


12) 


.25 2.5 


.720 1 
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OATA 



r i 



NAVAL P3ST GKAOUATC SCHCOl 



FILE : S\V 

li) .0625. 1.0*20 ,. lb 



13) 


• Uo ?:> 1 « 


► 0*2 0 


. 16 


14 ) 


.0 l f . 021/, .0235 






14) 


.0 l 


021 7 


.0235 


It) ) 


• Jd « • 1 2 5 » . J d 






15) 


• Od 


. 125 


.08 


L6P 


0 • , 4 . t !'• 04 






16) 


0. 


4. 


1.04 


17) 


0. f 4. , 1.04 






17) 


0. . 


4. 


1.04 


18) 


0 • f 50 0 • ,265. 






18) 


0. 


500 . 


23 5. 


19) 


0. t 500. ,342 • 






19 ) 


0. 


500. 


342. 


20) 


$ BLOCK G 






21 ) 


It / 






21 ) 


1 


7 




22 ) 


2 1 3 






22 ) 


2 


8 




23) 


3 1 9 






23) 


3 


9 




24 ) 


4,10 






24) 


4 


10 




25 ) 


5,11 






25) 


5 


11 




26) 


6, 12 






26) 


6 


12 




27) 


7,13 






27) 


7 


13 




28) 


8,14 






23) 


8 


14 




29) 


9,15 






29) 


9 


15 




30) 


$ 6 LOCK H 






31) 


11 






31 ) 


11 






32 ) 


24 






32) 


24" 






33) 


1 0.,2.5 






33) 


1 . 


0. 


2. 5 


34 ) 


25 






34) 


25 






35) 


•01 3,0. ,.13 






35) 


.0 13 


0. 


. 18 


36) 


26,27 






36) 


26 


27 




37) 


-1.0*20, 0., 0. 






37) 


— 1.0*20 


0. 


0 . 


38) 


1 9 






38) 


19 






39) 


- 1.0*20,0. , . 0722 






39) 


-1.0*20 


0. 


.0722 


40) 


20 






40) 


20 






41 ) 


— 1.0*20,0., • l*“t 







127 



FILE 



: SAV DATA Tl NAVAL PO ST.JR AOUAT E SCHOOL 



41) 

42 ) 

42) 

43 ) 


-1 .02 0 
30 

30 

• 7 f 0 • t 1.0*20 


0. 


.14 


43) 

44) 


. 7 

38, 38 


0. 


1.0*20 


44)« 
45 ) 


38 

0.0fJ.0,l.3 ,0. 


38 




45 ) 
*♦6 ) 


J. 0 

34,33 


0.0 


1.3 


46 ) 
47) 


3 4 

"“1.0,0. ,1. 


35 




47 ) 
46) 

48) 

49) 


-1. J 
36 

36 

2. ,0. , 1.0*7 


0. 


l . 


49 ) 
50) 
50) 
51 ) 


2. 

37 

37 

-1.02 0,0. ,0. 


0. 


1 . 0** 7 


51 ) 

52 ) 


-1 . C*2 0 
33 ,33 ,2 


0. 


0. 


52) 

53) 


33 

1. , . 1 , 1.0*20 


33 


2 


53) 

54) 

T IT LE 


1. 

EMD 


. 1 


1 .0*20 


AIK 


‘heated design - 


OPT I MIZE 


VOLUME 



CONTROL PARAMETERS; 

CALCULATION CONTROL, NCALC = 2 

NUMB 6r! OF \>L 03 AL DESIGN VARIABLES, NOV » 9 

INPUT INFORM AT IUN PRINT CDUE, IPNPUT = 0 

DEBUG PRINT CODE, I P DBS = 0 



CALCULATION CONTROL » NCALC 
VALUE MEANING 

1 SINGLE ANALYSIS 

2 OPT I MI Z AT I 0 R 

3 SENSITIVITY 

A TRU-VAk UBLE FUNCTION SPACE 

5 OPTIMU! SENSITIVITY 

6 APPROXIMATE OPTIMIZATION 



* * OPTIMIZATION INFORMATION 

GLOBAL VARIABLE iMUMOER OF OBJECTIVE = 16 

MULTIPLIER (NEGATIVE INDICATES MINIMIZATION) = -O.IOOCE 01 
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FILE 



: $AV 


OA TA 


ri 


NAVAL POST bR ADUA TE 


SCHOOL 


CCNRIN 


PARAMETERS (IF ZE-tO, C ON Ml N JEFAULT WILL 


OVER-R IJE) 




I PR IN 1 


I TM AX 


I CNJIR 


NS CAL ITRM LINJbJ 


NACMXl NFOG 


0 


ao 


0 


10 0 0 


1 1 


0 


FDCH 




FDCHM 


CT 


CTMIN 




0*0 




0.0 


0.0 


0.0 




CTL 




CTLMIN 


THETA 


PHI 




0.0 




0. J 


0.0 


0.0 




DEL FIN 




DAS F UN 


ALPH AX 


ABOBJi 




C.O 




0 .J 


0. 0 


0 . 0 




OES IGN 


VARI AOLE 


I MF JR HAT ION 







NON-ZERO HiriAL V4LJE A ILL OVER-RIDE MflDULE INPUT 



D. V. 


L U /* E * 


UP PE 4 INITIAL 








NU. 


8UJNJ 


SOUND VALUE 


SCALE 






1 


0.23200? uO 


J.2J2506 01 0.6 766JE 


00 0.0 






2 


0 • 2 b OOOE 00 


0.2 600 OG 01 0. 720IJE 


JO 0.0 






3 . 


0.623U0E- Ul 


0.110 JOE 16 0.160 JOE 


00 0. ) 






4 


0. L 00 JOc-0 1 


0.21 70 Jc-0 1 0.2 JbOJc 


-01 0.0 






5 


0.80 JOOE- 01 


0.125006 00 0.800 Jot 


-01 0.0 






6 


0.0 




0.4)00 06 01 0. 10400c 


Oi 0.0 






7 


0.0 




0.400006 01 0 . 104 ODE 


01 0.0 






8 


0.0 




0.500006 03 0. 235JJE 


03 0 .0 






9 


0.0 




0. 50000E 03 0.3420JE 


03 0.0 






DESIGN V ARI AdL 65 












D. V. 


GLCdA L 


MULTI PL7PIG 








ID 


NO. 


V4K. NU. 


FAC TUR 








1 


1 


7 


0. 100 OJE 01 








2 


2 


8 


0.100C0E 01 








3 


3 


9 ' 


0.100 COE 01 








4 


4 


10 


o.iooooe oi 








5 


5 


1 1 


0 *100006 01 • 








6 


6 


12 


O.IOOOOE 01 








7 


7 


13 


0.100 JOE 01 








8 


3 


14 


O.IOOOOE 0 1 








8 


9 


15 


0.100 JOE 01 








CCNSTRAINT INFORMATION 










THERE 


ARE 1 1 


CONSTRAINT 


SET S 










GLG8AL 


GLUbAL LINEAR LOWER NORM AL I Z AT ION 


U PPER 




ID 


VAR . 1 


VAR . 2 


ID 80 UN D 


FAC TOR 


SOUND 




1 


2<t 


0 


0 O.IOOOOE 01 


0. 10000? 01 


0 .2 5 00 ) d 


01 


3 


25 


0 


0 0. 1 30 000- 01 


0. 1 39J0E-0 L 


0 . 1 do j Jc 


00 


5 


26 


27 


0 -0. 1 1 000 E 16 


O.ilOUOE Lo 


0.0 




7 


19 


0 


0 -0. 1100-JE 16 


J . 1 1 0 J J E I o 


0 . 72200; 


-0) 


8 


20 


0 


0 -0. 1 10006 16 


0.1 100 06 16 


0. 140, JOE 


00 


9 


30 


0 


. 0 0. 70000E 00 


0. 70 000 c 3 j 


0.1 IOjOG 


lc 


10 


3d 


38 


0 0.3 


O.IOOOOE 00 


0. 13000 E 


01 
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DA T4 



Tl 



NAVAL POSTGRADUATE SCHOOL 



FILE: SAV 



12 


34 


35 


0 


-0. 100 00 G 


01 


0. IOOOOG 
0. 200 )0r. 


01 


0 . 100005 


01 


16 • 


36 


0 


0 


0. 2 00 001; 


01 


01 


0.10 )0 Jt 


06 


13 


3 7 


0 


0 


-0. 1 1000*5 


16 


0. 1 LUJOli 


16 


0 .J 




19 


33 


33 


2 


0. 1 00 001: 


01 


0* I OUOOt: 


00 


0. 1 iOJGE 


u 



TCTAL NUMBER OF CUN3TKAINEU PARAMETERS = 13 



* * ESTIMATED DAI A STCRAGE REQUIREMENTS 

REAL INTEGER 

INPUT EXECUTION AVAILABLE INPUT EX GLUT UN AVAILABLE 

106 650 600J • 76 146 1000 

1 AIR-COOLED HEAT EXCHANGER OPTIMIZATION 



INPUT DATA 



TUBESIDE PARAMETERS 

HASS FLOW R A TE= 133333.3 125 L3M/NR 

INLET TcUPEKATURE = 200.0000 DEG F 

OUTLET TEMPERATURE* 125 .0000 DEG F 

SPECIFIC FEAT* 1 . 00 00 BTU/LbM-F 



AIRSIDE PARAMETERS 



MASS FLCW SATE= 
INLET TEMPERATURE 

OUTLET temperatu* 
SPECIFIC FhAT = 
INLET PRESSURE* 



1190476. 0000 LBM/HR 

95.0000 DEG F 
= 1 3 0.0 COO JEG r 

0.24G0 GTU/LHM-F 
14.0000 PSI 



TUBE GEOMETRY 

TUBE INSIDE DIA.= 0.3700 INCHES 

TUBE OUTSIDE OIA.= 1.080) INCHES 

FIN HEIGHT= 0 .4600 INCHES 

FIN THICKNESS* G.OLOO INCHES 

FIN SPACING, CENT ER-T J-CEMT ER= 0. 1110 INCHES 



TUBE ARRANGEMENT 
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FILE: SAV 



OATA 



n 



NAVAL PO ST GRADUATE SCHCOL 



TRANS 4FRSE PITCH* 2.1250 INCHFS 

LONGUUjINAl Plfv-H* 2.1253 INCHF.S 

BANK HE HIM* 46.0000 INCHES 

BANK WIJTH* ' 2dtt .OOOJ INCHES 



INTEGEK PARAMETERS 

TYPE CF CROSS FUIw ARRANGEMENT* 7 

TYPE CF FIN PROFILE* I 

NUMBER OF RODS* A 

NUMBER OF PASSES* 2 



MISCELLANEOUS VARIABLES 

GIVEN FEAT TRANSFER RATE* IOOOOOOJ . 0000 8TU/MR 

THERMAL CONDUCTIVITY OF TH3S MATERIAL* 220.0000 3TU/HR-FT-F 

THERMAL CONOJCTIVITY OF' FIN MATERIAL* LIB. 0000 BTU/HR- FT-F 

GAS CONSTANT* 53.3400 FT-LBF/LBM-R 



penalty parameters 

12. CCCCCOOOO 0.0 0.0 



OPTIMIZATION PARAMETERS 

0.0 0.0 0.0 0.0 



* CC = O.IOOOOE 02 CMULT = 0.2C000E 01 CCMAX = 0. 10000E 04 
KOUNT* I 



* * CGKMIN DEIECTS INITIAL XID.GT.VU8U) 

XII) = 0.2J>jc-01 VOBIU = 0.2I70E-01 

XI I) IS SET EQUAL TO VUB(I) FOR I = 4 

KCUNT* 1 OEL * 0.77434E 01 OBJ* 0.753l0t 02 OBJ! * 0.723 C9R 02 
-0.435 IOE 00-0 . 42 5 (, o c 00-0. 205 E-0 3- 0 . 90000C 00-0.23553E 01 

— 0 • 8 44 04 £— J 1 -0 . o3 3 2 } T.“ J 1 -0 . 69 f')c. _• — Ci -0 .4 143 46 00 -0 .of') 3 4E 01 
-0.47/CiF 00-0. 13 2 6 7L JI-3.37124E J0-9.14932E jl -0 . 5 04 7 J E 00 
-0. 1 0044 E 03-0.909936 00-3.16773E 06 J.77484E 00 



0.0 


0. 0 


0.0 


0. 0 


0.0 


0.0 


0.0 


0. 0 


0.0 


0.0 


0 .0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


u. u 


0.77434E 01 
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Fite: SAV 
KOlINT = 



2 



□ ATA 



T 1 



NAVAL P3 ST CR ADUATE SCHOOL 



KOJNT = 2 OE L * 0. ? 2 1 26 F 01 UCU= 0.71452C 02 0QJ1 = 0.729226 02 

-0.454096 d'J-J. 4lH.l*>t 90-U. 132 JvC-Ol-O. 399606 00-0.101205 01 
-0.l7ftMF-d5-J.J772 1L -0 2-0. 42‘> 2 l E- 01 - ') .4 l .3 oil- JO-J. 65290c 01 
-0. 497/76 dO-U .16)756 01-0.0^2616 CO- J. I 500.) 6 Ol-O.OJOOOE 90 



-0 • 1 L i 5 9 £ 


03-0. 9 ?993 L 


00 —0 • 1 b 4 * 0 3 1 


Gb-0.2212bE 


JO 




0,0 


U. 0 


0.0 


J .0 




0.0 


0.0 


0.0 


0 .0 


0.0 




0.0 


0.0 


0.0 


0.0 


0.0 




0.0 


0.0 


0.0 


0 .0 


0.55358E 


01 




KO UN T = 


3 











KOUNT - 3 DEL= C.7C331E 00 QOJ= 0.71379E 02 OBJl = 0.722436 02 

-0.A5477E 00-0 . 4 L 6 c 00- 0 . 3 26 62 E - Cl-0 . 3 <>A6 S 6 00-0.10623E 01 
-0. 26009 E-0 1-0. 1 5 0 73 c -01 -0.031 70E-J2-0 .4141 IE 00-0 .oj40oE 01 
—0.4 76 l>c 00— 0.16091 E Jl-C.Jvdd/b 00-0 . I 4 u 'i5 t Jl-v>. 501 556 00 
-0.11315E 03 — 0. 9>9 4tl 6 00-0.164936 06-0 .7 00 3 IE-01 



0.0 




0.0 0.0 


0.0 


0.0 




0.0 




0.0 0.0 


0. 0 


0.0 




0.0 




0.0 0.0 


0.0 


0 .0 




0.0 

KO UN T = 


4 


0.0 0.0 


0. 48320 E 


01 




KOJNT = 




4 DEL* O.L C330E OL 


03J = 0 .71 241 E 02 U3J 1 = 


0.71464E 02 


-0.4544SE 


JO 


-O.MS^Oc 00-0. L3JiJ*E- 


01-0. 3981 *6 


OJ-O. 13 5 3 5 6 


01 


-0. 87023E- 


- I* 


-0 .321 80 £-02 -0. 3 40*+ 3 E- 


02-3 . tl3 J2E 


00-0.65 / 3 3E 


JL 


— 0.4 94 36E 


oO 


-0. Lb l U)l: JL-0.3 3J99E 


CO-J. L 50J0E 


31- 0.500 0 IE 


00 


-0. 113C9E 


03 


-U.S99SML UU-9. 16396E 


06-0 . 51651E- 


*31 




0.0 




0 • 0 O.*0 


0*0 


0.0 




0.0 




0.0 0.0 


0.0 


0.0 




0.0 




0.0 0.0 


0.0 


0 .0 




0.0 

KO UN T = 


5 


0.0 0.0 


0. 3 7 990 E 


01 





KOUNT = 


5* OE L = 


0.1 C 04 3 E 01 


03J-= 0.7L294E 02 J3J1 = 


0.7 1459 £ 


-0.4 54 5 1 E 


00-0.4L820C 


00-0 • L 6-*> 54E- 


Gl-O. 89835E 


00-0. 18 5 38 E 


01 


— 0. 56634E- 


0 j— J . 33o 32 b 


-02-0.3497 1E- 


02- J .91392 E 


00-0.637 33E 


0 1 


-0.4 94 36 E 


00-0 . 16 1 LO L 


J 1- 0 • 3 83 99fc 


00-0. 1 5000 c 


J1-0.5J003E 


00 


-0. 1 13G9E 


03-0.99993 fc 


00-0. L6 39 7 E 


06-0 . 502 - CE- 


•01 




0.0 


J .0 


0.0 


0.0 


0.0 




0.0 


0.0 


0.0 


0. 0 


0.0 




0.0 


0 .0 


0.0 


0.0 


0.0 




0.0 


0.0 


0.0 


0. 2 7942 E 


01 




KO UNT= 


6 











KOUWT= 6 DEL = 0.9412SE 00 C3J = 0.71 338E 02 03J1 = 0.71448E 02 

-0.45455E 0 0— t) .413136 00-0. 1 2 7 23E- 9 1- 0 . 6 59 73 E 00-0.11599 6 3 1 
-O.L 70 17 6-0 2-0 . 3 70 3 ■» 6-02-0 .33 0 4 t 2-02-0 . 4 13 7 JE 00-0 . 0 6 73 OE 01 
—0 .494366 OO-J.lollJc JI-0.3690JE CU-O.14099E 01-J.5J010E OO 
-0. 1 1309b 03-0.745936 00-0.16J69E Oo-O . 4 7364E- 01 



0.0 


0.0 


0.0 


0.0 


0 .0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 .0 


0.0 


0. 0 


0.0 


0.18529E 01 





132 



FILE: SAV 



DATA 



u 



NAVAL POSTGRADUATE SCHOOL 



KOUNT* 7 

KOUNT * 7 061 = 0. 1*33136 01 UHJ = 0.714026 02 D3J1 = 0, 

-0. 454566 00-0. 4 13 18 C JO-J. I 13 856- Cl-O. 8s896E 0J-0.l8o03E Ul 
-0 .2 123 71- 02- 0. J 3 7 o3l-) 2- ) . 39 l 20E- 22-0 I 1916 OO-O.o j 7 J3E 0 1 
-0.494 366 JJ-O.IoUjl Ol-C.JolOUt 00-0 . 140 >76 01-0.300136 00 
-0.113006 03-0.77998L 00-0.164006 06- 0 . 4 5 7956-0 1 



71445E 02 



0,0 


0.0 


0.0 




0.0 0.0 




0,0 


0.0 


0.0 




0.0 0.0 




0,0 


0. 0 


0.0 




0.0 0.0 




0,0 


O.J 


0.0 




0.210666-01 




KOUNT = 


d 










KOUNT- 


3 OE L- 


0. 1 75 17E 


01 


OBJ - 0.714336 02 OBJ l 


= 0. 


-0.45452E 


00-0. 4lH19t 


JO-O.l 1310E- 


91-0. 8 9.18 7 5 00-0.136066 


31 


-0.243<:4E- 


C2-0.4 I V .?£> E 


-J 2-0 .3^2 


2 -JE- 


02-0.413146 00-9.657336 


01 


-0.4943oE 


JO-O. 1 6 1 LOt 


01-9.389 OOF 


•JO-O. 14J99F: 01-0.500186 


00 


-0.1 131 JL 


23 -0. 9 >9 * j L 


JO-O .16 + 


03 E 


06-9.437335-0 1 




0.0 


0.0 


0.0 




0.0 0.0 




0.0 


0.0 


0.0 




0.0 0.0 




0.0 


0. 0 


0.0 




0.0 0.0 




0.0 


0.0 


0.0 




-0. 173066 01 





KOUNT= 



KOUNT= 0 06 L= 0.169336 01 03J= 0.71555E 02 UBJi = 0.714466 02 

-0.454496 00-0.410206 JO-9. 1 12 oo5- 01-0. 6 93.3 7E 00-9.166036 01 
-0.26 4 596-02-0. 43 044 6- 0 2-0 . 39 8 05C- 02-0 . 4 1 3 946 00-0 .65 7 3 36 01 
-0.494366 00-0.161101- 01-0.33 0096 00-0. 149936 Jl-0. 500166 00 
— 0 . I 1 3 1 2E 03-0.99993 6 00-0 . 1 6+ C56 Jo-3.424586-0 1 



0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 .0 


0.0 . 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

KO UN T = 


0.0 

10 


0.0 


-0. 34239 E 01 





KOUNT - 
-0.452726 
-0. 140 77E 
-0.49460 6 
-0.113446 
0.0 
0.0 
0.0 
0.0 

KO UN T = 



10. 06 L= 0.254596 01 U3J = 0.71 5306 02 OBJ 1 = 0.715996 02 
00-0.413016 JO- 0. 135246-01-0. b 986 5 6 JO-O. 136896 01 
-01-0.13301L-91-0.65 7646-02-0.414116 00-0.o>7026 01 

00-0. 1499 25 01-9.300846 00 
06 0.318246-01 

0.0 0.0 0.0 0.0 

0.0 0.0 0 . 0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 -0.883026 00 

11 



■J 0-0 . I <j 1 0 86 91-0.38)246 
03-0. 997' 7 8c OO-J.lo-3 226 



KOIJN T = ll 06 L= 0.257546 01 OSJ= 0 .716096 02 0QJ1 = 0.7 1596E 02 

-0.45271 6 90-0.413916 00-0.131 056-01-0. 8936 9= 00-0.166396 01 
-0. 14149E-01-0. 162 /36-JI-0.56 1996-02-0 .4141 1 6 JJ- 0.65701 E 01 
—0.49481 6 0 0-0 . 16 10 5 c 01-0.38 0 2 r ’ 6 OJ-O. 149926 01-0.6 )0846 JO 
— 0 .11 3446 C3-0.999 6 )0-J.lo522E 06 0.321925-01 



0.0 


0 .0 


0.0 


0.0 


0 .0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.16924E 01 
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FILE: SAV 
KOUN T= 



DATA 



TI 



NAVAL POSTGRADUATE SCHOOL 



12 

KOUNT= 12 OEl = '3.259366 01 
-0.45271 I JJ-J.41843( JJ-O. 12 Jo 76 
-0. 141 34C-J 1-0. l6 2o 7 t- ) l-).5o J 436 
— 0 • 49 4 6 1 1 JO-O.lolO/l: 31-0.334256 
-0.1134*6 J3-0.4 444Jt 30-0 .103 22 E 



0 .0 


J. 0 


0 .0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 .0 


0.0 

KOUN T= 13 


0.0 


0.0 


KCU NT = 


13 OEL = 


0.171316 Oi 



-0.45173 6 JO- J . 4 l 92 Oh JJ-0.2 0r>34 6 
-0. I 7735 u-0 1-0. I 3 53 4 6- Jl- 3.1 15 826 
-0.49t346 30-0. l»l Opc 01 - 0 . 3 39V4 6 
-0.113726 OJ-U./'r^d: )0- 0 . 1 64 9 /C 



0.0 


0. 0 


0 .0 


0.0 


0.0 


0.0 


0 .0 


0.0 


0 .0 


0.0 


0.0 


0.0 


KOUN I s 


14 





K CU NT - 14 OFL= O.lo9536 01 

—0.451 736 00-0. 414296 00-0.20^636 
-0.1 77SC6-J1-0. 13 504 E-0 1-0. 11582C 
-0.49484E OO-O.lolO^C Ol-O.ih 344 fc 
-0.110721; 00-0.49 9'4de 00-0.i6t97E 



0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


KGUNT= 


15 




K CJ NT = 


15 • OEL= 


0.16739E 01 


-0*431 7 6 E 


0O-Q.41929C 


00-0.2 C3 1 7 E 


-0.1 7835E 


-Oi-J. 13 504;; 


-0 1-0. 1 1 5 32 E 


-0.49-UWE 


00-0 . 1 *j 1 Oj t 


31 -0.33040 E 


-0.1 1372E 


03-0.999906 


00- 0.164 9 76 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 .0 


0.0 


KOUNT- 


16 




KQJ NT = 


16 061= 


0.3 32^56 01 


-0.431736 


00-0.41 929F 


JQ-0 .2 u ? 2 6 6 


-0.1 73GCE 


-01 - J. l 


-J 1-0 .113 326- 


-0.<t9434E 


00-0. 16105b 


31-0.33 9 4*6 


-0.11371b 


03-9.9^99^ l 


JO- 0 • 1 G^ 9 76 


0.0 


0. 0 


0.0 


0.0 


0.0 


0.0 


0.0 


J. J 


0.(3 


0.0 


0 .0 


0.0 



UO J - 0.71 691 £ 02 08 J 1 = 0.715 946 02 

Jl-O. 39 4/11 00-0 . I 16 Od E 01 
02 -0 .41411c 00-3 .6-3 7006 01 
00-0 . I *♦ 992 E 01 -0.6 JO rt 36 00 
06 0.3 242'JE-Ol 



0.0 0.0 

0.42860E 01 

OB J = 0.71 7 39E 32 30 J1 = 0.71 776E C 2 

01-3.397)41: 00-0. 1 36056 01 
Cl-3.4143 76 00-0. 65671E 01 
00-3 . I 49396 01-0. 501 06E 30 
C6-0. I 0 73 5S-01 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.25679E 01 



UBJ = 0.717576 02 DQJ 1 = 0. 71 775E C2 
01-3 .397486 00-0 . 1 36366 01 
01-0.4 1-0 76 03-0 . & 56 7 l 6 31 
30-3.149596 01 - 0.301 Oo 6 JO 
06-0. 105956-01 

0.0 0.0 

3.0 • 0.0 

0.0 0.0 

0.8726CE 00 



OBJ = 0.71774E 02 ' J B J 1 = 0.71 7746 C2 

01-3 . B 9 7 J 7 6 J 3-0 . I 36 3c> 6 01 
01-0.41^076 GO-0. 696716 01 
00-3 .144896 01-0.301076 00 
06-0. 104626-01 

0.0 0.0 

3.0 0.0 

0.0 0.0 

-3.801296 00 



OB J = 0.717996 32 0BJ1 = 0.717736 C2 
31-3.6974 86 JO-0 .1 56 38 6 01 
01-0.414J76 00-0.656716 01 
00-3.1 44896 01-0.531 076 00 
06-0. I 0405 6-01 

0.0 0.0 

0.0 0.0 

0.0 0.0 

-0.413036 01 
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FILE: SAV 



OATA 



T 1 



NAVAL PQST3RA0UATE SCHOOL 



KOUN T = 



17 



KOUNT = 1 7 061= J.2124JE 01 010= 0.718036 02 OBJl 

-0.45134E 00-0 . 4 1 Vt 7 6 JO- ) . 22 2 l 3L - 01 -0. 867 hi fc 0 J-J . lJo >2 E 
-0. 10 7B2c-0 1-) .l^l Jc-Jl-O. 1 18 2 IE- ) 1- J. AIM 2 C 00- J. 6 os *06 
-0.4T4W21- 00-0. 1<> U-»t 1 1-0 . J;i-) 57E CO-O. lAOOHI: 01-J.3J1 18E 
-0.113311 OJ-O.O'.VMt JU-J. l0 5 2Vc 06 0.663 7<j 6-02 

0.0 0.0 0.0 0 .0 

0.0 0.0 0.0 O.J 

0.0 0.0 0.0 0.0 

0.0 0.0 -0.2 000 7E 01 



= 0.71924E 02 
Jl 
01 

00 



o.o 
0 .0 
0.0 
0.0 

KOUN T •= 



10 



K OUNT = 13 OEL = 0.21353E 01 OBJ = 0.71317E 02 OUJi = 0.71823E 02 

-0.451J4E 00 -0.410. 72 00-0 . 22 l 32E- Cl -0. 5 9 77 0£ JJ-J.l0oV2E 01 
-0. 198 0 IE- J 1-0 . I oto 13 l - 0 1-0. 1 15 l >c-01-0. 4 14 12C 00-0 .086 60c Jl 
— 0 .4 94 02 E 0 )-') . lu l.)>L U— 0.3 3 5*3 76 0J-0.14OHF 01-0.30116E JO 
-0.113616 OJ-0.4J843L 00-0. lo5 20 E 06 0.66 7456-02 



0.0 


0 . 0 


0.0 


0.0 


0 .0 




0.0 


0.0 


0.0 


0.0 


0.0 




0.0 


0.0 


0.0 


0.0 


0 .0 




0.0 


0.0 


0.0 


0.12912E 


00 




KOUN T = 


19 










KQUNT - 


19 DE 1= 


3.42976F 01 


D BJ = 0.71 8 37E 02 03 J l = 


0.71822E 


•0.451 34E 


00—0 .41946c 


00-0. 220 79 E - 


Cl - 0. 8 9 7 7 9 E 


00-0 . 1 0 o 9 2 E 


91 


-0. 1 18 10 b 


-0 1-3. I6o02b 


- J 1 —0 • 115 1 9 E- 


31-0. 4 1 41 2 E 


JO- 0.65660b 


01 


-0.49*92E 


00-0. lo LJ**E 


J 1-0.38957E 


C J —0 . 1 4 9 8 8 E 


01-0 .50 USE 


00 


-0.11381c 


03- 0. 59993b 


30-0. 168 29 E 


06 0. 67151 E- 


•02 




0.0 


0.0 


0.0 


0.0 


0 .0 




0 .0 


0. 0 


0.0 


0.0 


0.0 




0.0 


0.0 


0.0 


0.0 


0 .0 




0.0 


0.0 


0.0 


0 .4426 8 E 


01 




KQ UN T = 


20 










KOUNT - 


20 OE L= 


0.43011b 01 


OOJ= 0.719 6*6 1 02 G3J1 


0.71821E 


-0.46134E 


00-0.41946 E 


0 0-0.220 3qE- 


01-0. 8 9 7 9 0 E 


00-0 . 186 9 2 E 


)1 


-0. 11603E- 


-0 1-0 . 1d5 7 >L 


-01-0. 115 19c- 


0 1 - v ) • 41412c 


00- 0 . 8 56 o 1 E 


01 


-0.494 92E 


00-0 . 16 10 tt 


) 1-0 . 3 69d /E 


00—0 • 1 4 93 n E 


01-0 .DO 1 1 BE 


00 


-0. 1 138 1 E 


03-0. S559JL 


10-0. 163 29 E 


Jo J.672U4E- 


02 




0.0 


0.0 


0.0 


0.0 


0 .0 




0 .0 


0. J 


0.0 


0.0 


0 .0 




0.0 


0.0 


0.0 


0. 1 


0 .0 




0.0 


0. 0 


0.0 


0 • 87273 E 


01 





KQUNT= 20 OEL= 0.43011E 01 C0J= 0. 71 3946 02 08J1 -= 

1 OPTIMIZATION 3 c SO LT S 



OBJECTIVE FUNCTION 

GLOBAL lQCATUN 16 FJNCTICN VALUE 0.71321E 02 



DESIGN VARIABLES 
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FILE: SAV DATA Ti 



NAVAL POSTGRADUATE SCMCOL 



10 

1 

2 

3 

4 
6 
6 

7 

8 
9 




GLCIbAL 
VAR. NU. 

7 

8 
9 

10 

ii 

L2 

13 

14 

15 



LOxER 

30 JMO 

0 • 2320 jE 03 
0.2 50 JOE 00 
0. o25 006-0 L 
0. 999996*02 
0. 79999E-01 
O.G 
0. 0 
0.0 
0.0 



VALUE 

0.63H306 00 
0.5/5596 00 
0. 1298 9b OC 
u. 2 9 7 / >6-0 1 
0. 79 >9 >6-01 
0 • d3 7 36 6 00 
0. 1022 56 OL 
0. 1 3 60 7 E 03 
0.236036 03 



UPPER 
00 UN C 

0.2 325 Cf* 01 
0 . 250006 01 
o.iioor- u 
3.21 7006-01 
0. 12 50 06 00 
0.39999c 01 

0. -0000c 91 
0.4 999 5F 03 
0.499996 03 



DESIGN CONSTRAINTS 



1 



GLOBAL 



ID 


VAR. 


1 


24 


3 


25 


5 


26 


6 


27 


7 


19 


8 


20 


9 


30 


10 


38 


12 


34 


14 


35 


16 


36 


13 


37 


19 


33 



LQw c R 

. 3UJN0 

0.1 00 CO 6 01 
0. ldOJJE-Oi 
-0. 1100 JE 16 
-0.110006 16 
-0. 1 1 000 6 lb 
-0.110006 16 
0 • 700 OJ 6 00 
0.0 

-o. uoooe o i 
- 0 . 1000)6 01 
0 . 200 00 6 0 1 
-0.110006 16 
0.10U006 01 



VALUE 

0.14513E 01 
0.1 83976-01 
-0.18692b 90 
-0.193036-02 
0 . 7 100 36- j i 
0.130596 00 
0.909806 00 

0. o566 06 )0 
0.610436 00 
0.4yd62c )0 
0 . 2 2 96 ^6 03 
-0.165296 05 
0. 99933 6 00 



UPPER 

BOUND 

0.2 5 0 COE 01 
0. UOOOE 00 
0.0 
0.0 

0. 72 2 C 06-0 1 
0. 140 COE 00 
o.iiooor ib 

0.1 30 )06 01 
0. 100006 01 
O.lOOCOr: 01 
0. 1 00006 03 
0.0 

o. uoooe 16 



OUTPUT DATA 



TUBESIDE PARAMETERS 

MASS FLOW RATE= 133333.3125 L3M/HR 
INLET TEMPERATURE = 200.0000 06 G F 

OUTLET T6 MP6 R4 TUR c = 125.D0J0 JEG F 

SPECIFIC FcAT - 1.0000 6TU/L8M-F 



A1RSI0E PARAMETERS 



136 



FILE: SAV 



oata 



ri 



NAVAL POSTGRADUATE SCHOOL 



MASS FLOW RATE* 
INLCT TEMPERATURE* 
OUTLET TEMPERATURE* 
SPECIFIC HE A [ = • 
INLET PRESSURE* 



1190476.0000 LOM/HR 

95. COCO OtG F 

130.0 000 OEG F 
0. 2 A 00 OTU/LOM-F 

14.0000 PS1 



TUBE GEJMETRY 

TUBE INSIDE DIA.* 0.5383 INCHES 

TUBE UUfSlOi: 0 1 A. = 0.5756 INCHES 

FIN HEIGHT* 0.1299 INCHES 

FIN THICKNESS* 0.020? INCHES 

FIN SPACING, CENTER-TU-CCNTE.R* 0.0800 INCHES 



. TUBE ARRANGEMENT 



TRANSVERSE PITCH* ■ 1.0223 INCHES 

LONGITUDINAL PITCH* 0. 8374 INCHES 

BANK HEIGHT* 236.0935 INCHES 

BANK rtlOTH* 136.0704 INCHES 



THETA* 0.6556 

DELPHI* -0.3920 



INTEGER PARAMETERS 

TYPE CF CECSS FLO* ARRANGEMENT* 7 

TYPE CF FIN PROFILE* l 

NUMBER OF RUNS* 4 

NUMBER OF PASSES* 2 

VERTICAL ROWS CONTAIN 229 TUBES 



MISCELLANEOUS VARIABLES 

GIVEN HEAT TRANSFER RATE* 10000000.0000 B TU/HR 

THERMAL CCNuJC T IV IT i CF TUJt MATERIAL* 220. 0000 BTU/H.R-FT-F 

THERMAL CONDUCTIVITY OF FIN MATERIAL* 113. 0000 BTU/HR-F T-F 

GAS CCNSTANT* 53.340c FT-LBF/L3M-R 



PENALTY PARAMETERS 

12.COOOOOOOO 0.0 0.0 
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FILE: SAV 



DATA 



NAVAL PJSTGRAO’JATE SCHOOL 



T I 



OPTIMIZATION PARAMETERS 

0.0 0. 0 



0.0 



0.0 



OTHER VALUES 

HEAT TRANSFER RATtS 
ODO T 1 = 99 S 32 EG. JOOO QTU/HR 

Q0GT2* 999S993. JOOJ 

OCOT 3= v959<; ^.oouo 

HEAT TRANSFER DIFFERENCES 
OE LCI = 672 J.OOOJ 3 TU / H^ 

DELU2= -o 71 3. UOUO 

DEL Ca= 4.UJU0 



PENAL 1Y = 



11.9919 



CONSTRAINTS 
DIAMETER RATIO* 

TUBF THICKNESS* 

TDUCHN* - 0.1 867 

T CUCH L= - 0.0020 

PROFH - - 70.0000 F 

PROFC= - 30.0000 F 



VALUES TO BE OPTIMIZED 
VOLUME* 71.32 LA CJ.FT. 

A RE A= 10105.4219 oO.FT. 

A IR HP - 92 . 7 29 o HP 

AIR PRESSJRE DROP* 

TUBESIDE PkESSJRE DROP* 

OBJ ECT IVE FUNC T I ON* 

PROGRAM CALLS TO ANALIZ 



1.4513 

0.0134 INCHES 



0.0710 PSI 

0.1334 PSI 

71.8214 



ICALC CALLS 

1 i 

2 1150 

3 1 



138 



APPENDIX D 



HEDSUP PROGRAM LISTING 
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PITCHN - DISTANCE BETWEEN CENTERS OF ADLACENT TUBES NORMAL TO FLOW 
(TRANSVERSE PITCH), INCHES 



OOOOOOOOOOOOOOOOOOOCOOOOOCOOOOOOOOOCCOOOOOOOOOCO 
nO h- ao cr o <n m in so h cc O' o *-» cm ro < m vo r- co o> o p-i <n m <■ in so h oo o o ^ eg m <■ m so h oo O' o — « n g co 
^ vf <• -4" inm m mm m in m m m no o o -o vC o o >o Hr- r- h r * h h h r-r*^ cc x cc oc ao cc coco cn oo cr O' O' o 
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WRITE <6,130) AMDOT, ARTMPl , AR1MP2 , AGP , PRE S 1 
WRITE (6,140) TUBEID»TUBE03,FINHT»FINTH,FINSP 
WRITE (6,150) PI TCHN,PITCHL ,BANKH,BANKW 
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